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1.  Introduction 

Increased  expression  of  eIF4E  has  frequently  been  reported  in  breast  cancers  and  is  thought 
to  make  fundamental  contributions  to  disease  development  and  progression1.  Increased 
eIF4E  activity  acts  to  enhance  the  translation  of  cancer-related  transcripts  that  contain 
highly  structured  5’  untranslated  regions  (UTRs)  in  their  mRNAs.  Over-expression  of 
eIF4E  has  been  shown  to  correlate  with  poor  prognosis  in  breast  cancer2  therefore  the  level 
of  eIF4E  over-expression  has  been  studied  extensively  as  a  prognostic  marker  with  some 
success.  Furthermore,  eIF4E  is  an  established  target  for  cancer  therapy3  and  clinical  trials 
of  the  efficacy  and  safety  of  cancer  therapeutics  that  target  eIF4E  have  been  carried  out, 
again,  with  some  success.  However,  it  is  clear  that  eIF4E  expression  does  not  equate  to 
eIF4E  activity,  since  eIF4E  activity  is  additionally  regulated  by  a  family  of  binding 
proteins,  the  4E-BPs,  that  bind  to  and  inhibit  eIF4E  activity4.  Our  hypothesis  was  that 
analysis  of  eIF4E  activity  in  individual  breast  tumours,  as  opposed  to  eIF4E  expression, 
gives  improved  prognostic,  predictive  and  biological  understanding  of  individual  breast 
cancers,  and  overall  insights  into  the  mechanisms  of  breast  carcinogenesis. 

In  order  to  test  this  hypothesis  we  aimed  to  investigate  eIF4E  activity  in  breast 
carcinogenesis  using  a  novel  cross-disciplinary  approach.  First,  we  aimed  to  determine 
expression  levels  of  eIF4E  and  its  regulators  in  breast  tumours,  and  to  estimate  the 
individual  contributions  of  each  factor  to  eIF4E  activity  by  relating  each  expression  level  to 
cancer  survival.  This  approach  would  allow  us  to  build,  in  collaboration  with 
mathematicians,  an  equation  relating  the  levels  of  these  factors  together  to  estimate  eIF4E 
activity.  Next  we  aimed  to  validate  these  mathematical  relationships  in  tissue  culture 
systems  where  we  would  be  able  to  measure  both  expression  levels  and  eIF4E  activity. 
Finally,  we  aimed  to  determine  whether  measures  of  eIF4E  activity  in  cell  lines,  and 
estimates  of  eIF4E  activity  in  tumours,  would  act  as  predictive  markers  for  cellular 
responses  to  eIF4E  directed  therapies. 
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2.  Body 

2.1  Examination  of  archival  breast  tumours  (Statement  of  Work  sections  Al,  Bl,  B4). 

We  have  determined  expression  levels  of  eIF4E,  4E-BP1,  4E-BP2  and  phospho-4E-BPl 
Thr37/46  (p-4E-BPl)  by  immuno-histochemistry  (IHC)  in  breast  tumours  from  a  cohort  of 
424  patients,  supported  by  extensive  clinical  background  and  follow  up.  We  have 
established  the  relationship  between  expression  of  these  markers  and  tumour  grade,  size 
and  type.  Moreover,  we  have  combined  the  insights  gained  from  each  of  these  markers  into 
an  estimate  of  eIF4E  activity.  We  have  developed  a  mathematical  function  that  relates 
expression  levels  of  eIF4E,  4E-BP1,  4E-BP2  and  p-4E-BPl  in  breast  tumours  to  survival. 
This  function,  V,  provides  additional  prognostic  insights  when  compared  to  examination 
of  eIF4E  expression  levels  alone.  This  variable  can  be  described  as  X-B1/4+PB1/2-B2/4, 
where  X,  Bl,  PB1  and  B2  represent  eIF4E,  4E-BP1,  p-4E-BPl  and  4E-BP2  levels 
respectively.  This  work  has  been  published5  (manuscript  attached  within  appendices). 

2.2  Determination  of  expression  levels  of  eIF4E  and  its  regulators,  and  estimation  of 
eIF4E  activity,  in  cell  lines  (SoW A3,  B2). 

Next,  we  were  interested  to  determine  the  validity  of  our  estimate  of  eIF4E  activity.  We 
aimed  to  study  tissue  culture  cells  in  which  we  would  be  able  to  determine  both  expression 
levels  of  eIF4E  and  its  regulators  (this  section),  and  actual  eIF4E  activity  (section  2.3). 

We  determined  relative  expression  levels  of  eIF4E,  4E-BP1,  4E-BP2  and  p-4E-BPl  in  a 
panel  of  human  cell  lines  (SW480,  Caco-2,  U2020,  H1299,  A549,  HB2,  MDA-MB-231, 
MCF-10A,  MCF7)  using  Western  blotting  and  densitometry.  Two  alternative  cell  lysis  and 
normalisation  methods  were  used:  RIPA  lysis  normalized  to  total  protein  (Figure  1),  and 
Laemmli  lysis  normalized  to  expression  of  beta-actin  (data  not  shown).  In  order  to 
reproduce  the  quantification  that  was  used  for  breast  tumours  (section  2.1),  expression  of 
the  markers  in  the  cells  lines  was  also  determined  using  IHC.  We  found  no  correlation 
between  assessments  of  relative  expression  by  Western  and  IHC  (Figure  1).  It  is  worth 
noting  that  assessment  by  IHC  in  cell  lines  proved  particularly  problematic.  We  found 
marked  heterogeneity  of  expression  throughout  the  cell  populations.  In  addition,  cellular 
morphology  was  poorly  maintained  making  differentiation  between  cytoplasmic  and 
nuclear  staining  difficult. 
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Figure  1  Expression  levels  of  eIF4E  and  its  regulators  as  determined  by  Western  blot  and 
densitometry  (filled  bars)  or  IHC  (open  bars)  do  not  correlate.  Expression  levels  of  eIF4E, 
4E-BP1,  p-4E-BPl  and  4E-BP2  in  a  panel  of  human  cell  lines  were  determined  by 
densitometric  analyses  of  Western  blots  of  RIP  A  extracts  normalized  to  total  protein,  or  by 
semi-quantitative  scoring  of  intensity  and  the  proportion  of  cells  showing  positivity  after 
immuno-histochemical  staining  of  agarose  encased  and  formalin  fixed  paraffin  embedded 
cells.  Data  were  normalised  to  values  from  SW480  cells  to  allow  comparisons. 

Next,  we  determined  estimates  of  eIF4E  activity  in  these  cell  lines  by  calculating  z  (see 
section  2.1)  using  the  expression  data  from  Western  blots  or  IHC.  This  is  not  a  trivial 
procedure  for  the  data  from  Western  blots.  Densitometry  values  from  Western  blots  were 
first  mapped  to  values  that  were  within  an  appropriate  range  for  use  with  the  z  function  -  ie 
from  0  to  7.  However,  simply  linearly  transforming  densitometry  values  to  fit  within  this 
range  defining  the  lowest  as  0  and  the  maximum  value  as  7  does  not  represent  a  likely 
solution  with  only  nine  cell  lines,  since  scores  of  0  and  7  occur  rarely  in  tumours.  We 
transformed  the  maximum  value  in  the  cell  lines  to  the  most  common  positive  value  in  the 
tumours  (eIF4E:  5;  4E-BP1:  6;  p4E-BPl:  5;  4E-BP2:  5)  and  have  maintained  the  fold 
differences  to  this  for  the  other  cell  lines.  As  expected,  considering  the  differences  in  the 
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assessments  of  expression  levels  by  Western  and  IHC,  there  was  no  significant  correlation 
between  the  z  values  determined  from  these  two  datasets  (Figure  2). 


Figure  2  Estimated  eIF4E  activities  based  on  data  from  Western  blots  (filled  bars)  or  IHC 
(open  bars)  do  not  correlate.  eIF4E  activities  were  estimated  from  expression  data  for 
eIF4E,  4E-BP1,  p-4E-BPl  and  4E-BP2  using  the  function  z5. 

2.3  Determination  of  eIF4E  activity  in  cell  lines  (SoW A4,  A5,  B3). 

The  cellular  role  of  eIF4E  is  to  bind  to  mRNA  caps  allowing  recruitment  of  eIF4F,  and 
subsequently  the  translational  machinery.  The  complex  formed  scans  linearly  along  the 
5’UTR  until  an  initiation  codon  in  good  context  is  encountered,  at  which  point  further 
elements  of  translational  machinery  are  recruited  and  protein  synthesis  starts.  Some 
5’UTRs,  frequently  those  on  cancer-related  transcripts,  contain  a  high  degree  of  secondary 
structure  that  can  inhibit  cap-recognition  and  translational  scanning.  Increased  eIF4E 
activity  is  thought  to  reduce  the  effects  of  5’UTR  structure  by  enhancing  cap-recognition 
and  scanning,  therefore  increasing  translation  of  these  specific  oncogenic  transcripts1. 
Having  produced  two  different  estimates  of  eIF4E  activities  in  cell  lines  (Figure  2),  we  next 
wished  to  establish  the  accuracy  of  the  estimates.  Therefore,  we  aimed  to  measure  actual 
eIF4E  activity  experimentally.  Two  techniques  to  measure  eIF4E  activity  were  used. 


Firstly,  we  used  a  fluorescent  reporter  assay.  We  have  previously  shown  that  a  5’UTR 
expressed  from  the  human  AXIN2  gene  contains  a  sixty  nucleotide  sequence  that  is 
predicted  to  form  a  stable  stem-loop  structure6.  This  structure  fits  the  criteria  associated 
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with  UTRs  that  determine  differential  translational  efficiencies  in  response  to  changes  in 

n 

eIF4E  activity  ,  while  lacking  other  motifs  known  to  allow  regulation  of  translation  (e.g. 
upstream  AUG  codons  or  binding  sites  for  trans-acting  proteins).  We  now  wished  to 
examine  whether  the  translational  efficiency  defined  by  this  UTR  would  respond  to 
experimentally  induced  changes  in  eIF4E  activity,  and  could  therefore  be  used  to  measure 
endogenous  cellular  eIF4E  activity.  The  sequence  was  cloned  upstream  of  the  GFP  reading 
frame  in  an  expression  vector.  MCF7  cells  were  transiently  transfected  with  an  equal  copy 
number  of  vectors  to  allow  expression  of  GFP  mRNAs  with  either  a  control  non-regulatory 
5 ’UTR  or  this  structured  5 ’UTR,  along  with  either  empty  expression  plasmids  or  plasmids 
allowing  eIF4E  over-expression.  GFP  protein  expression  was  measured  by  flow-cytometry 
and  GFP  mRNA  expression  was  measured  by  qPCR  allowing  determination  of  relative 
translational  efficiencies  for  each  GFP  message  (Figure  3A).  Western  blot  analyses  were 
used  to  confirm  expression  of  exogenous  eIF4E  in  the  appropriate  co-transfected  cells 
(Figure  3B).  The  translational  efficiency  of  the  control  reporter  was  not  significantly  altered 
by  eIF4E  over-expression  (compare  lanes  1  and  2),  demonstrating  that  eIF4E  over¬ 
expression  did  not  cause  a  general  enhancement  of  translation.  As  previously  reported6,  the 
structured  5’UTR  conferred  repression  of  translation  (compare  lanes  1  and  3;  p=0.002). 
Critically,  this  repression  was  overcome  by  exogenous  eIF4E  (compare  lanes  3  and  4; 
p=0.002),  resulting  in  translation  with  the  same  efficiency  as  messages  lacking  inhibitory 
5 ’UTRs.  We  concluded  that  this  reporter  did  indeed  respond  to  changes  in  eIF4E  activity. 
Relative  translational  efficiencies  specified  by  this  eIF4E-responsive  5’UTR  were 
determined  in  the  panel  of  cell  lines.  Cells  were  transiently  transfected  with  vectors  to 
allow  expression  of  GFP  mRNAs  with  control  or  the  structured  5’UTR  as  before,  and 
translational  efficiencies  were  determined  (Figure  4A).  A  range  of  translational  efficiencies 
was  seen,  with  A549  cells  determining  the  lowest,  and  HB2  cells  the  highest  (Figure  4B). 

Note  that  we  have  now  employed  the  translational  efficiency  assay  developed  here  for 
further  published  studies  of  UTR  function8  (manuscript  attached  within  appendices). 
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1  2  3  4  -  -  +  elF4E 

elF4E:  +  + 
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Figure  3  Translational  efficiency  specified  by  a  structured  5’UTR  reporter  responds  to 
changes  in  eIF4E  activity.  A)  Reporters  were  constructed  to  express  mRNAs  containing  the 
GFP  reading  frame  preceded  by  a  control  5’UTR  lacking  regulatory  motifs  (control)  or  a 
sequence  predicted  to  form  a  stable  stem-loop  structure  (structured).  MCF7  cells  were 
transiently  transfected  with  equal  copy  numbers  of  either  control  or  structured  reporters 
along  with  either  empty  expression  vector  (-)  or  vector  to  allow  over-expression  of  eIF4E 
(+).  GFP  protein  and  mRNA  were  quantified  by  flow-cytometry  and  real-time  PCR 
respectively.  Translational  efficiency  (protein  synthesised  per  unit  mRNA)  is  presented 
relative  to  the  control.  Data  points  represent  means  (+/-standard  deviations)  of  technical 
triplicates  within  a  representative  experiment.  B)  Expression  of  exogenous  proteins  was 
confirmed  by  Western  blot  analysis  using  anti-HA  (exogenous  GFP  and  eIF4E  both  include 
the  HA  epitope). 


A  SW480  Caco2  U2020  B 


Figure  4  Translational  efficiencies  specified  by  a  structured  5’UTR  reporter  vary  in 
different  cell  lines.  A)  Cells  were  transiently  transfected  with  equal  copy  numbers  of 
plasmids  to  allow  expression  of  transcripts  with  the  GFP  reading  frame  preceded  by  either  a 
control  5’UTR  lacking  regulatory  motifs  (con)  or  a  sequence  predicted  to  form  a  stable 
stem-loop  structure  (struc).  GFP  protein  and  mRNA  were  quantified  by  flow-cytometry  and 
real-time  PCR  respectively.  Translational  efficiency  (protein  synthesised  per  unit  mRNA) 
is  presented  relative  to  the  control.  Data  points  represent  means  (+/-standard  deviations)  of 
technical  triplicates  within  a  representative  experiment.  A  minimum  of  two  independent 
experiments  were  performed.  B)  Translational  efficiencies  of  transcripts  with  structured 
5’UTRs  in  the  9  cells  lines. 
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Secondly,  an  assay  that  allows  quantification  of  eIF4E  that  is  capable  of  binding  to  the 
mRNA  cap  was  used.  Methyl-7-GTP-sepharose  was  used  to  purify  cap-binding  eIF4E  from 
lysates  of  each  cell  line  containing  equal  masses  of  total  protein;  eIF4E  recovered  was 
quantified  using  Western  blotting  and  densitometry  (Figure  5A).  These  values  represent 
relative  quantification  of  eIF4E  that  can  bind  mRNA  caps  within  these  cell  lines,  which 
could  be  thought  of  as  ‘eIF4E  activity’.  However,  we  found  that  different  amounts  of  4E- 
BP1  co-purified  with  eIF4E  in  each  cell  line  (Figure  5B).  Consequently,  we  concluded  that 
differing  proportions  of  the  eIF4E  that  was  active  in  terms  of  cap-binding  in  each  cell  line 
would  not  have  been  able  to  assemble  eIF4F  and  initiate  translation.  Therefore  we  have  not 
used  this  assay  as  a  physiologically  relevant  measure  of  eIF4E  activity  in  terms  of  the 
oncogenic  function  of  enhancing  translation  of  specific  transcripts. 


Figure  5  Cap-binding  assays  determine  eIF4E’s  cap-recognition  activity  but,  since  4E-BP1 
is  co-purified,  this  may  not  relate  to  eIF4E’s  ability  to  assemble  active  translation 
complexes.  Proteins  were  purified  on  methyl-7-GTP-sepharose  from  cell  lysates  containing 
equal  masses  of  total  protein.  Purified  eIF4E  (A)  and  4E-BP1  (B)  were  quantified  by 
Western  blot  and  densitometry. 


Estimates  of  eIF4E  activity  that  were  based  on  expression  levels  of  eIF4E  and  its  regulators 
(section  2.2),  and  the  measures  of  eIF4E  activity  from  the  translational  reporter  assay  (this 
section)  were  compared  (Figure  6).  No  significant  correlations  were  observed  between 
estimates  (black  and  open  bars)  and  the  actual  measurement  of  eIF4E  activity  (grey  bars) 
(Spearman’s  rank  correlations:  Western  derived  z  with  measured  activity  r=0.367,  p=0.33; 
IHC  derived  z  with  measured  activity  r=-0.161,  p=0.678).  The  lack  of  correlation  with  the 
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estimate  derived  from  IHC  expression  data  may  reflect  the  inaccuracy  of  IHC  assessment  in 
the  context  of  cell  lines.  The  estimate  derived  from  Western  blot  data,  in  fact,  appears  to  be 
surprisingly  accurate  in  five  of  the  nine  cell  lines;  this  may  suggest  that  further  regulatory 
factors  not  taken  into  account  in  the  current  Western  analyses  may  have  critical  roles  in  the 
other  four  cell  lines. 
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Figure  6  Estimated  and  experimentally  determined  eIF4E  activities  do  not  correlate.  eIF4E 
activities  were  estimated  using  the  z  function  (see  Figure  2)  on  data  derived  from  Western 
blots  (black  bars)  or  IHC  (open  bars).  eIF4E  activity  was  experimentally  determined  by 
analyses  of  the  translational  efficiency  specified  by  a  structured  5’UTR  (grey  bars;  see 
Figure  4B).  Data  were  normalised  to  values  from  SW480  cells  to  allow  comparisons. 


2.4  Analyses  of  the  therapy  predictive  value  of  estimates  and  measures  of  eIF4E  activity 
in  cell  lines  and  in  tumours  (SoW Al,  B4). 

Rapamycin,  and  derivatives  such  as  RAD001  (Everolimus),  have  been  tested  as  cancer 
therapeutics  with  some  success9.  The  drugs  inhibit  the  activity  of  the  protein  kinase 
mammalian  Target  Of  Rapamycin  (mTOR).  It  is  well  established  that  mTOR  activity  is  up- 
regulated  in  many  cancers10.  Critical  mTOR  targets  are  4E-BP1  and  2,  which  -  when 
hyperphosphorylated  by  mTOR  -  are  unable  to  interact  with  and  inhibit  eIF4E.  Thus  cells 
in  which  mTOR  activity  is  high  are  thought  to  have  high  eIF4E  activity.  Our  hypothesis 
was  that  high  eIF4E  activity  would  be  associated  with  high  sensitivity  to  mTOR  inhibitors, 
as  the  cells  would  require  continued  mTOR  activity  to  maintain  their  eIF4E-dependent 
aberrant  growth  and  survival. 
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To  test  this  hypothesis  in  cell  lines,  we  treated  our  panel  of  cell  lines  with  rapamycin  and 
determined  their  relative  sensitivities.  Cells  were  treated  with  different  doses  of  rapamycin 
and  proliferation/survival  was  determined  relative  to  control  treated  cells  using  MTT  assays 
(Figure  7  A).  Sensitivities  to  the  highest  dose  are  shown  in  Figure  7B.  A  range  of 
sensitivities  was  seen,  with  a  three-fold  difference  between  the  most  sensitive  (MCF7)  and 
most  resistant  (MDA-MB-231). 


Figure  7  Cell  lines  have  different  sensitivities  to  rapamycin.  A)  Cells  were  treated  with 
either  control,  or  different  doses  of  rapamycin  and  growth/proliferation  was  monitored 
using  MTT  assays.  MTT  readings  after  48  hours  are  shown  relative  to  control.  B)  Relative 
sensitivities  to  rapamycin  are  shown;  these  are  the  %  reductions  in  growth/proliferation 
caused  by  lOOnM  rapamycin.  Data  points  represent  means  (+/-  standard  deviations)  from 
five  technical  replicates  within  a  representative  experiment.  A  minimum  of  two 
independent  experiments  were  performed. 

Sensitivities  were  compared  with  the  eIF4E  activities  determined  using  the  reporter  assay 
(section  2.3).  Initially,  we  analysed  this  relationship  using  Spearman’s  rank  correlation 
coefficient;  we  found  a  strong  and  significant  positive  association  (r=0.72;  p=0.037). 
However,  the  correlation  was  particularly  evident  in  8  of  the  9  cell  lines;  if  MCF7  cells, 
which  were  more  sensitive  to  rapamycin  than  predicted,  were  excluded  from  the  analysis 
the  strength  and  significance  of  the  relationship  was  increased  (r=0.83;  p=0.015).  Similarly 
in  linear  regression,  a  highly  significant  relationship  was  seen  when  MCF7  cells  were 
excluded  from  the  analysis  (Figure  8;  p=0.0037).  We  concluded  that  eIF4E  activities  as 
determined  by  our  reporter  assay  did  predict  sensitivity  to  rapamycin  in  8  of  the  9  cell  lines. 
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Figure  8  Experimentally  determined  eIF4E  activities  correlate  with  sensitivities  to 
rapamycin.  Data  from  Figure  4B  (x-axis)  and  Figure  7B  (y-axis)  were  plotted  for  each  cell 
line.  Linear  regression  was  performed  to  determine  the  relationship  in  the  cell  lines 
excluding  the  outlier,  MCF7  cells;  the  linear  model  is  shown  as  a  line  (black)  with  95% 
confidence  intervals  (grey  lines)  (p=0.0037). 

We  have  also  examined  whether  our  estimate  of  eIF4E  activity,  the  z  function,  provides 
predictive  insights  with  respect  to  response  to  eIF4E-targetted  therapy  in  clinical  breast 
cancers.  Tumour  samples  from  23  breast  cancer  patients  treated  with  at  least  10  days  neo¬ 
adjuvant  RAD001  were  obtained  both  pre-treatment  (core  biopsies  taken  for  diagnosis)  and 
post-treatment  (surgical  resection  samples)  from  collaborators  in  Edinburgh  (Prof  John 
Bartlett,  Edinburgh  Cancer  Research  UK  Centre).  Samples  were  stained  by  IHC  for  eIF4E, 
4E-BP1,  4E-BP2  and  p-4E-BPl  as  before.  Scoring  of  expression  levels  remains  underway. 
Two  independent  observers  will  score  the  samples  and  then  consensus  scores  will  be 
reached  for  those  where  there  is  a  discrepancy.  Currently,  only  scores  from  one  observer 
are  available.  Firstly,  it  is  clear  that  RAD001  significantly  reduced  p-4E-BPl  levels  (in 
21/23  cases).  Therefore,  it  appears  that  RAD001  acted  as  expected,  to  inhibit  mTOR 
activity  and  therefore  reduce  mTOR-dependent  phosphorylation  of  4E-BP1.  However,  the 
situation  is  more  complex,  since  total  4E-BP1  was  also  reduced  in  10  of  these  cases  - 
therefore  whether  the  pool  of  hypo-phosphorylated  4E-BP1  that  was  free  to  inhibit  eIF4E 
was  increased  is  unclear.  Furthermore,  changes  in  levels  of  eIF4E  or  4E-BP2  were  also 
commonplace.  In  order  to  estimate  the  changes  in  eIF4E  activity  induced  by  RAD001 
treatment  the  values  of  z  were  calculated  both  before  and  after  treatment.  The  z  value 
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decreased  in  only  16/23  cases.  Our  collaborators  have  determined  the  responses  of  the 
patients  to  the  RAD001  treatment  in  terms  of  changes  in  the  proliferation  indexes  of  their 
tumours  (assessed  by  Ki67  staining).  When  our  assessment  of  eIF4E  activities  is  complete, 
we  will  analyse  whether  estimated  eIF4E  activity  at  diagnosis,  or  the  change  in  eIF4E 
activity  after  RAD001  treatment,  correlates  with  response  to  the  drug;  thus  whether  z  acts 
as  a  predictive  or  pharmacodynamic  marker. 
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3.  Key  Research  Accomplishments 

1)  We  have  established  that  mathematically  combining  assessments  of  expressions  of 
eIF4E-regulators  with  assessments  of  expression  of  eIF4E  in  clinical  tumours  provides 
improved  prognostic  insights  over  examination  of  eIF4E  alone,  presumably  as  this 
combinatorial  function  estimates  eIF4E  activity.  This  work  has  been  published5. 

2)  We  have  demonstrated  that  an  experimentally  determined  measure  of  eIF4E  activity 
predicts  cellular  sensitivity  to  rapamycin  (manuscript  under  revision),  suggesting  that 
measurement  of  eIF4E  activity  within  tumours  could  be  used  as  a  predictive  marker  for 
therapy  with  RAD001  or  other  rapamycin-derivatives  in  cancer. 

3)  We  are  in  the  process  of  testing  whether  our  combinatorial  estimate  of  eIF4E  activity 
tumours  can  be  used  as  a  predictive  marker  for  neo-adjuvant  RAD001  therapy  in  breast 


cancer. 
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4.  Reportable  outcomes 

This  work  has  led  directly  to  two  published  papers  (attached  in  appendices)  and  one 
manuscript  under  revision: 

Coleman  LJ,  Peter  MB,  Teall  TJ,  Brannan  RA,  Hanby  AM,  Honarpisheh  H,  Shaaban  AM, 
Smith  L,  Speirs  V,  Verghese  ET,  McElwaine  JN,  Hughes  TA  (2009)  Combined  analysis  of 
eIF4E  and  4E-binding  protein  expression  predicts  breast  cancer  survival  and  estimates 
eIF4E  activity.  Br  J  Cancer  100,  1393-9 

Smith  L,  Coleman  LJ,  Cummings  M,  Satheesha  S,  Shaw  SO,  Speirs  V,  Hughes  TA  (2010) 
Expression  of  estrogen  receptor  beta  isoforms  is  regulated  by  transcriptional  and  post- 
transcriptional  mechanisms.  Biochem  J  429,  283-290 

Satheesha  S,  Coleman  LJ,  Ingram  N,  Madhok  B,  McElwaine  JN,  Hughes  TA. 

Translational  efficiencies  defined  by  a  structured  5’  untranslated  region  predict  cellular 
sensitivity  to  rapamycin.  Under  revision 
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5.  Conclusions 

We  have  tested  and  supported  our  main  hypothesis:  analysis  of  eIF4E  activity  in  breast 
tumours  gives  improved  prognostic  insights  into  breast  cancer  over  analysis  of  only  eIF4E 
expression. 

We  have  also  demonstrated  that  measurement  of  eIF4E  activity  provides  a  predictive 
marker  for  response  of  cell  lines  to  treatment  with  the  drug  rapamycin.  Work  is  still 
underway  to  establish  whether  our  clinically  practicable  method  of  estimating  eIF4E 
activity  in  tumour  samples  provides  a  predictive  marker  for  response  of  breast  tumours  to 
neo-adjuvant  treatment  with  the  rapamycin-derivative  RAD001. 

Finally,  we  have  developed  a  generally-applicable  cross-disciplinary  methodology  for 
combining  the  prognostic/predictive  insights  from  individual  biomarkers  into  estimates  of 
the  activities  of  biological  pathways.  This  approach  could  be  applied  to  estimate  pathway 
activities  and  therefore  predict  response  for  any  molecularly-directed  cancer  therapy.  The 
most  obvious  example  is  the  her2  pathway  that  is  targeted  by  herceptin  (trastuzumab).  This 
therapy  is  assigned  on  the  basis  of  over-expression  of  the  biomarker  and  molecular  target 
her2.  However,  a  significant  proportion  of  patients  fail  to  respond  to  the  therapy,  while  it 
has  been  reported  that  some  her2-negative  tumours  do  respond.  It  is  likely  that  an  estimate 
of  her2  activity,  based  on  combined  assessments  of  expression  of  her2  and  its  direct 
regulators,  would  provide  an  improved  predictive  biomarker. 
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Increased  eukaryotic  translation  initiation  factor  4E  (elF4E)  expression  occurs  in  many  cancers,  and  makes  fundamental  contributions 
to  carcinogenesis  by  stimulating  the  expression  of  cancer-related  genes  at  post-transcriptional  levels.  This  key  role  is  highlighted  by  the 
facts  that  elF4E  levels  can  predict  prognosis,  and  that  elF4E  is  an  established  therapeutic  target.  However,  elF4E  activity  is  a  complex 
function  of  expression  levels  and  phosphorylation  statuses  of  elF4E  and  elF4E-binding  proteins  (4E-BPs).  Our  hypothesis  was  that  the 
combined  analyses  of  these  pathway  components  would  allow  insights  into  elF4E  activity  and  its  influence  on  cancer.  We  have 
determined  expression  levels  of  elF4E,  4E-BPI ,  4E-BP2  and  phosphorylated  4E-BPI  within  424  breast  tumours,  and  have  carried  out 
analyses  to  combine  these  and  relate  the  product  to  patient  survival,  in  order  to  estimate  elF4E  activity.  We  show  that  this  analysis 
gives  greater  prognostic  insights  than  that  of  elF4E  alone.  We  show  that  elF4E  and  4E-BP  expression  are  positively  associated,  and 
that  4E-BP2  has  a  stronger  influence  on  cancer  behaviour  than  4E-BPI.  Finally,  we  examine  elF4E,  estimated  elF4E  activity,  and 
phosphorylated  4E-BPI  as  potential  predictive  biomarkers  for  elF4E-targeted  therapies,  and  show  that  each  determines  selection  of 
different  patient  groups.  We  conclude  that  elF4E’s  influence  on  cancer  survival  is  modulated  substantially  by  4E-BPs,  and  that 
combined  pathway  analyses  can  estimate  functional  elF4E. 
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The  eukaryotic  translation  initiation  factor  4E  (eIF4E)  has  key 
roles  in  carcinogenesis  (De  Benedetti  and  Graff,  2004).  eIF4E  is 
often  overexpressed  in  carcinoma  cells  as  compared  to  equivalent 
normal  epithelium  in  many  tumour  types  including  breast 
(Kerekatte  et  al ,  1995),  lung  (Rosenwald  et  al ,  2001)  and  colon 
(Rosenwald  et  al ,  1999).  The  oncogenic  role  of  this  overexpression 
has  been  shown  by  various  experimental  observations;  for 
example,  forced  eIF4E  overexpression  within  many  cell  types 
leads  to  transformation  (De  Benedetti  and  Graff,  2004),  and  within 
transgenic  mice  increases  incidence  of  multiple  tumour  types 
(Ruggero  et  al>  2004).  eIF4E  has  at  least  two  normal  cellular 
functions.  First,  it  is  an  essential  component  of  the  multimeric 
factor  eIF4F,  which  initiates  cap-dependent  translation  -  the 
mechanism  responsible  for  most  protein  synthesis  (Gray  and 
Wickens,  1998).  eIF4E’s  role  is  to  bind  to  mRNA  caps  allowing 
recruitment  of  eIF4F,  and  subsequently  the  translational  machin¬ 
ery.  The  complex  formed  scans  linearly  along  the  5'  untranslated 
region  (UTR)  until  an  initiation  codon  in  good  context  is 
encountered,  at  which  point  further  elements  of  translational 
machinery  are  recruited  and  protein  synthesis  starts.  Second, 
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eIF4E  regulates  expression  of  some  genes  by  controlling  nuclear 
export  of  their  transcripts  (Culjkovic  et  al>  2007),  a  function  that 
also  requires  eIF4E’s  cap-binding  activity  (Culjkovic  et  al ,  2005). 
Under  most  normal  conditions,  availability  of  active  eIF4E  is 
thought  to  be  rate  limiting  for  both  functions.  One  might  expect 
a  general  translational  stimulation  to  result  from  the  increased 
eIF4E  expression  in  cancers,  on  account  of  enhanced  mRNA  cap 
recognition,  yet  effects  of  eIF4E  overexpression  are  more  subtle. 
Approximately  10%  of  mammalian  transcripts  have  5'UTRs  that 
may  form  complex  secondary  structures  that  reduce  the  abilities  of 
both  eIF4F  to  bind  to  mRNAs  and  the  translational  machinery  to 
scan  5'UTRs  (Pesole  et  al ,  2001);  the  result  is  that  these  transcripts 
are  translated  inefficiently  (Hughes,  2007).  The  majority  of  human 
transcripts  with  these  inhibitory  5'UTRs  code  for  growth  or 
cancer- associated  proteins  (Kozak,  1991).  Increased  eIF4E  is 
thought  to  reduce  the  effects  of  5'UTR  structure  by  enhancing 
cap-recognition  and  scanning,  therefore  increasing  translation  of 
these  specific  oncogenic  transcripts  (De  Benedetti  and  Graff,  2004). 
Similarly,  increased  eIF4E  expression  enhances  nuclear  export  of  a 
set  of  transcripts  associated  with  oncogenesis  (Culjkovic  et  al , 
2007).  As  a  consequence  of  this  central  role,  eIF4E  is  an  established 
target  for  cancer  therapy  (Smolewski,  2006;  Graff  et  al ,  2008). 

The  importance  of  eIF4E  in  cancer  has  been  underlined  by  the 
fact  that  eIF4E  expression  levels  can  be  used  to  determine 
prognosis.  Cases  in  which  eIF4E  is  highly  overexpressed  tend  to 
have  poor  prognoses  (Li  et  al ,  1998).  A  substantial  confounding 
factor  is  that  eIF4E  expression  does  not  equate  to  eIF4E  activity, 
thereby  making  interpretation  of  potential  influences  of  eIF4E 
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levels  difficult.  eIF4E  activity  is  a  complex  function  of  eIF4E 
expression  and  expressions  and  activities  of  eIF4E-binding 
proteins  (4E-BP1,  2  and  3)  that  bind  to  and  inhibit  eIF4E  (Richter 
and  Sonenberg,  2005)  (Supplementary  Figure  SI).  Activity  is 
further  regulated  by  phosphorylation  of  4E-BP1  (and  other  4E-BPs 
it  is  assumed),  with  only  hypophosphorylated  forms  being  able 
to  inhibit  eIF4E.  Additional  regulation  occurs  by  differential 
phosphorylation  of  eIF4E  itself,  although  there  are  conflicting 
reports  as  to  how  this  influences  activity  (Scheper  and  Proud, 
2002).  The  result  is  that  high  expression  of  eIF4E  may  not  lead  to 
high  eIF4E  activity  if,  for  example,  hypophosphorylated  4E-BP1 
were  also  highly  expressed.  Many  cancer- related  signalling  path¬ 
ways,  including  PI3K  and  p38,  converge  to  regulate  eIF4E  and 
4E-BP  phosphorylation;  therefore,  eIF4E  activity  seems  to  be  a  key 
cancer-signalling  node  (Polunovsky  and  Bitterman,  2006).  Here, 
we  have  tested  the  hypothesis  that  combined  analyses  of 
expressions  and  phosphorylation  states  of  eIF4E,  and  its  regulators 
allows  greater  understanding  of  eIF4E  activity  and  its  influence  on 
cancer  than  examination  of  eIF4E  expression  alone. 

MATERIALS  AND  METHODS 
Patients 

Ethical  approval  was  obtained  (Leeds  East  05/Q1206/136).  Archival 
cancer  tissue  and  data  were  obtained  for  424  patients  diagnosed  at 
LTH  NHS  Trust  from  1983-2006.  Tissue  microarrays  (TMAs)  were 
constructed  containing  0.6  mm  cores  selected  from  representative 
tumour  areas  as  determined  by  a  consultant  breast  histopatho- 
logist  (AMS)  from  H&E  stained  sections.  Survival  periods  - 
overall:  initial  diagnosis  to  death;  disease-free:  initial  diagnosis  to 
the  diagnosis  of  recurrence/metastasis;  disease-specific:  initial 
diagnosis  to  death  after  recurrence  or  metastasis  (cancer-specific 
death  confirmed  in  most  cases). 

Westerns  and  immunohistochemistry 

MCF7  and  MDA-MB-231  cells  were  cultured/transfected  as  earlier 
(Johnson  et  al ,  2008;  Maraqa  et  al ,  2008).  An  eIF4E  expression 
vector  was  obtained  from  John  Blenis  (Harvard  Medical  School) 
and  Nahum  Sonenberg  (McGill).  Western  analyses  were  carried 
out  as  earlier  (Maraqa  et  al ,  2008)  using  the  reagents  in 
Supplementary  Table  SI.  TMA  sections  of  5  /urn  were  dewaxed 
and  blocked  in  hydrogen  peroxidase  block  (20  min).  Antigens  were 
retrieved  and  stained  as  described  in  Supplementary  Table  SI,  and 
as  used  elsewhere  (Zhou  et  al ,  2004;  Dutton  et  al ,  2005;  Lee  et  al , 
2005;  Engelman  et  al>  2008).  Envision  detection  was  used  (DAKO, 
Glostrup,  Denmark).  Negative  controls  (primary  antibodies 
omitted)  were  included  in  each  immunohistochemistry  (IHC) 
batch;  in  addition,  adjacent  normal  epithelium,  lymphocytes  and 
blood  vessel  endothelium  served  as  internal  controls.  Controls 
were  performed  for  p4E-BPl  antibodies  in  which  sections  were 
pretreated  with  Lambda  Phosphatase  (Nebraska,  NE,  USA).  Cores 
were  scored  for  immunoreactivity  by  two  or  more  individuals  (LJC, 
TJT,  ETV  and  RAB),  taking  into  account  the  average  intensity  and 
percentage  of  positively  stained  tumour  cells  (as  used  earlier  for 
eIF4E  (Zhou  et  al ,  2006)).  Staining  intensity  scores  (0  no  staining, 
1  weak,  2  moderate  and  3  strong)  were  added  to  percentages 
positively  stained  scores  (1  <5%;  2,  6-25%,  3,  26-75%  and  4 
>  75%),  giving  totals  of  0  or  2  -  7.  Consensus  scores  were  determined 
for  cores  with  different  initial  scores,  and  all  scoring  was  overseen  by 
a  consultant  breast  histopathologist  (AMH). 

Mathematical  analyses 

Data  were  analysed  using  Kaplan -Meier  survival  curves.  Depen¬ 
dence  on  prognostic  indicators  was  determined  using  Cox  propor¬ 
tional  hazards  models;  significance  values  relate  to  likelihood  ratio 


tests  of  the  null  hypothesis  that  indicators  do  not  effect  hazard 
rates  (Cox  and  Oakes,  1984).  SPSS  (SPSS,  Chicago,  IL,  USA)  and 
the  Statistics  Toolbox  in  MATLAB  (ecdf.m  and  coxphfit.m) 
(MathWorks,  Natick,  MA,  USA)  were  used.  Tests  were  two  sided 
and  P<0.05  was  considered  significant. 

RESULTS 

Antibody  validation 

We  have  used  IHC  to  determine  expression  levels  in  breast 
tumours  of  the  main  regulatory  molecules  of  the  eIF4E  pathway 
-  namely,  eIF4E,  4E-BP1,  4E-BP2  and  phosphorylated  4E-BP1 
(Thr37/46)  (termed  p4E-BPl).  We  have  not  examined  4E-BP3 
because  it  is  not  thought  to  have  a  role  in  breast  (Poulin  et  al , 
1998),  or  phosphorylated  forms  of  eIF4E,  as  their  influences  on 
the  activity  and  in  cancer  remain  uncertain  (Scheper  and  Proud, 
2002;  Salehi  and  Mashayekhi,  2006;  Buxade  et  al ,  2008).  First, 
we  optimised  the  antibody  use  on  archival  breast  tissue.  We 
established  that  antibodies  were  specific  for  their  antigens  using 
western  blots  against  lysates  of  breast  cancer  cell  lines  (Supple¬ 
mentary  Figure  S2A).  In  addition,  we  showed  phospho-specificity 
of  antibodies  against  p4E-BPl  by  carrying  out  IHC  on  serial  tissue 
sections  with  and  without  pretreatment  with  protein  phosphatase 
(Supplementary  Figure  S2B). 

Patient  cohort  and  immunohistochemistry 

Tissue  micro-arrays  containing  samples  from  424  breast  tumours 
were  established,  supported  by  detailed  clinicopathological  data 
(Supplementary  Table  S2).  The  cohort  included  a  wide  range  of 
patient  and  tumour  characteristics,  with  mean  patient  follow-up  of 
91.9  months.  We  carried  out  IHC  for  eIF4E,  4E-BP1,  4E-BP2  and 
p4E-BPl  on  TMA  sections  and  assessed  immunoreactivity  within 
tumour  cells,  taking  into  account  the  proportions  of  cells  staining 
positively  and  average  intensity,  giving  scores  of  0  (negative)  or 
2-7  (positive).  Representative  staining  patterns  are  shown 
(Figure  1,  and  at  higher  magnification  in  Supplementary  Figure 
S3).  Tumour  stroma  and  normal  tissue  were  negative  for  eIF4E, 
4E-BP1  and  4E-BP2,  whereas  very  occasional  low-intensity  staining 
for  p4E-BPl  was  noted  in  normal  epithelial  cells.  Staining  was 
generally  cytoplasmic,  although  nuclear  staining  was  noted  in  a 
minority  of  cases  (Supplementary  Figure  S4);  this  was  separately 
analysed  and  was  found  not  to  be  of  prognostic  value  and  is  not 
discussed.  As  expected,  data  were  not  available  for  some  patients 
because  of  the  TMA  core  loss  during  processing,  a  well-recognised 
occurrence,  therefore,  data  for  all  four  antigens  were  available  for 
only  282  patients.  The  full  range  of  scores  were  observed  for  each 
antigen  (Figure  II).  It  was  notable  that  staining  was  most 
frequently  not  detectable  for  p4E-BPl.  Others  have  reported  more 
frequent  expression  of  p4E-BPl  (Zhou  et  al ,  2004),  therefore  we 
carried  out  IHC  for  an  alternative  p4E-BPl  species  (Ser65);  we 
found  immunoreactivity  with  this  antibody  to  be  similarly 
infrequent  (see  discussion). 

Expressions  of  eIF4E  and  4E-BPs  correlate  with  grade 

Associations  between  antigen  expressions  and  a  wide  range  of 
clinicopathological  parameters  were  examined.  No  correlations 
were  found  with  nodal  status,  tumour  size  or  histological  type. 
Weak  positive/borderline  no  correlations  were  found  with  oestro¬ 
gen  receptor  a  status  and  eIF4E  expression  (Spearman’s 
p  coefficient  0.21;  P<  0.001)  and  4E-BP2  (0.22;  P<  0.001),  but 
not  with  4E-BP1  or  p4E-BPl.  Strong  correlations  between 
expression  of  markers  and  tumour  grade  were  found.  eIF4E 
expression  (split  into  three  classes,  0-3,  4-5  and  6-7)  was 
positively  associated  with  grade  (y2-test,  P  =  0.011),  whereas 
expression  of  both  4E-BPs  was  negatively  associated  with  grade 
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Figure  I  The  full  range  of  expression  intensities  and  proportions  for  eukaryotic  translation  initiation  factor  4E  (elF4E),  elF4E-binding  protein  (4E-BP)  I , 
4E-BP2  and  p4E-BP  I  occur  within  breast  cancers.  (A-H)  Representative  tumour  tissue  microarray  cores  showing  immunoreactivity  as  labelled.  These  cores 
were  scored  A  4,  B  7,  C  3,  D  7,  E  3,  F  7,  G  3  and  H  7.  (I)  Histograms  showing  distributions  of  immunohistochemistry  scores  within  breast  cancers.  Scores 
(x-axis)  and  numbers  of  cores  assigned  to  each  score  (y-axis)  are  shown.  Percentages  of  the  cohort  are  given  above  the  bar  for  each  score. 


(4E-BP1  P  =  0.002;  4E-BP2  P  =  0.029).  p4E-BPl  was  positively 
associated  with  grade  (P  =  0.012).  A  positive  association  between 
eIF4E  expression  and  grade  has  been  reported  earlier  (Li  et  al , 
2002).  Correlations  for  other  markers  were  consistent  with  their 
influences  on  carcinogenesis  being  through  the  eIF4E  pathway;  4E- 
BPs,  eIF4E  inhibitors,  were  negatively  associated  with  grade, 
whereas  4E-BP1  phosphorylation,  which  would  relieve  4E-BP1- 
induced  inhibition  of  eIF4E,  was  positively  associated. 

High  expression  of  eIF4E  correlates  with  poor  prognosis 

Kaplan -Meier  survival  analyses  were  used  to  determine  survival 
with  respect  to  eIF4E.  Analyses  were  carried  out  with  expression 
divided  into  IHC  scores,  although  scores  of  0,  2  and  3  were 
combined  as  each  individual  group  was  small,  for  overall  survival 
(OS),  disease-free  survival  (DFS)  and  disease-specific  survival 
(DSS)  (Figure  2A-C).  High  eIF4E  scores  were  indicative  of  poor 
prognosis.  Prognosis  seemed  to  worsen  with  each  increasing  eIF4E 
score  for  OS,  whereas  patterns  for  DFS  and  DSS  suggested  weaker, 
but  still  detectable,  influences  of  individual  scores  with  an  overall 
grouping  into  two  classes  (0-5  good  prognosis;  6  or  7  poor 
prognosis).  We  have  examined  relationships  between  eIF4E 
expression  and  survival  using  Cox  regressions.  We  have  included 
either  eIF4E  scores,  or  eIF4E  expression  dichotomised  arbitrarily 


or  as  suggested  by  the  apparent  bimodal  distribution  seen  above, 
and  modelled  these  with  respect  to  OS,  DFS  and  DSS.  Models  that 
most  accurately  reflected  the  data  included  eIF4E  scores  rather 
than  dichotomised  data,  showing  the  value  of  scoring  proportion 
and  intensity  of  positive  tumour  cells.  In  these  models  each 
increase  of  1  in  eIF4E  score  gave  increases  in  hazard  ratios 
(HRs)  of  1.22  {P  =  0.004),  1.3  (P  =  0.008)  and  1.33  (P  =  0.005)  for 
OS,  DFS  and  DSS,  respectively.  Thus  individuals  with  scores  of 
7  have  DFS  HRs  of  6.15  (95%  CIs:  3-12,  P  =  0.008)  as  compared 
with  individuals  with  scores  of  0.  The  prognostic  value  of 
eIF4E  has  been  reported  earlier  as  independent  of  grade/nodal 
status  in  breast  cancer  (Li  et  al ,  1998);  we  have  examined 
independence  from  the  Nottingham  Prognostic  Index  (NPI),  which 
takes  account  of  tumour  size,  grade  and  lymph  node  status 
(Haybittle  et  al ,  1982).  In  multivariate  Cox  regressions  the 
prognostic  value  of  eIF4E  was  independent  of  NPI  with  eIF4E 
remaining  significant  for  OS  (NPI  P<0.00001;  eIF4E  P  =  0.02), 
DFS  (NPI  P<  0.00001;  eIF4E  P  =  0.045)  and  DSS  (NPI  P<  0.00001; 
eIF4E  P  =  0.029). 

Prognostic  value  of  4E-BP1,  4E-BP2  and  p4E-BPl 

Kaplan -Meier  survival  analyses  were  also  used  to  determine 
survival  with  respect  to  the  other  markers.  We  present  data 


©  2009  Cancer  Research  UK 


British  Journal  of  Cancer  (2009)  1 00(9),  1393-1 399 


Clinical  Studies 


elF4E  activity  predicts  breast  cancer  survival 

LJ  Coleman  et  al 


1396 


Disease-specific  survival  (years) 


D  4E-BP1  E  4E-BP2 


F  phospho-4E-BP1 


Disease-free  survival  (years) 


Figure  2  Expression  of  eukaryotic  translation  initiation  factor  4E  (elF4E),  but  not  elF4E-binding  protein  (4E-BP)  1 , 4E-BP2  or  p4E-BPI ,  is  associated  with 
prognosis.  (A-C)  Kaplan -Meier  survival  analyses  for  overall  (A),  disease-free  (B)  and  disease-specific  survival  (C)  for  patient  groups  with  tumours  with 
differing  elF4E  scores.  (D-F)  Kaplan-Meier  survival  analyses  for  disease-free  suwival  for  patient  groups  with  tumours  with  differing  scores  for  4E-BPI  (D), 
4E-BP2  (E)  or  p4E-BPI  (F).  As  scores  of  0,  2  and  3  were  relatively  rare,  these  have  been  grouped  together.  Censoring  ticks  have  been  omitted  for  clarity. 


for  DFS  (Figure  2D-F),  and  for  OS  and  DSS  (Supplementary 
Figure  S5).  These  antigens  provided  little  prognostic  insight  and 
we  were  unable  to  construct  significant  Cox  equations  to  model 
their  individual  survival  influences.  In  the  case  of  4E-BP2,  there 
was  a  nonsignificant  trend  for  high  scores  to  associate  with  good 
prognosis. 

Mathematical  modelling  of  influences  of  4E-BPs:  4E-BPs 
modify  eIF4E  activity  and  provide  additional  prognostic 
insights 

Next  we  examined  influences  of  4E-BPs  in  the  context  of  eIF4E 
expression.  We  necessarily  restricted  these  analyses  to  the  282 
patients  for  whom  scores  of  all  four  antigens  were  available.  Using 
this  dataset,  models  including  solely  eIF4E  expression  gave  HRs  of 
1.21  (P  =  0.011),  1.24  (P  =  0.035)  1.27  (P  =  0.02)  for  OS,  DFS  and 
DSS,  respectively,  for  each  increase  of  1  in  eIF4E  score.  We  were 
unable  to  construct  significant  Cox  models  on  the  basis  of 
combinations  of  4E-BPs  without  including  eIF4E,  suggesting  that 
eIF4E  is  their  critical  effector.  We  found  models  combining 
expression  of  4E-BP1  or  p4E-BPl  (as  IHC  scores,  or  dichotomised 
into  two  groups)  with  eIF4E,  provided  little  additional  prognostic 
value  over  that  found  with  eIF4E  alone.  However,  including 
4E-BP2  in  a  model  for  OS  enhanced  the  model  significantly  with 
hazard  increasing  by  1.28  with  each  point  increase  in  eIF4E  score 
(P  =  0.005)  and  decreasing  0.11  with  each  increase  in  4E-BP2  score 
(P  =  0.02).  We  refined  this  by  combining  eIF4E  and  4E-BP2  scores 
into  a  single  non-linear  variable  in  which  high  levels  of  eIF4E  or 
4E-BP2  act  to  increase  or  decrease  the  value  respectively  (achieved 
using  max  (0,X-B2/3.5),  where  X  and  B2  represent  eIF4E  and  4E- 
BP2  scores).  This  variable,  termed  y,  predicted  survival  more 
accurately  than  examination  of  eIF4E  alone;  each  increase  of  1  in 
y  carried  HRs  of  1.32  (P=  0.0003),  1.32  (P  =  0.013)  and  1.36 
(P  =  0.006)  for  OS,  DFS  and  DSS,  respectively.  We  also  investigated 


Cox  models  including  expression  of  all  four  antigens.  In  order  to 
combine  terms,  we  considered  models  that  included  each  variable 
separately  and  found  their  relative  effect  on  HRs  using  maximum 
likelihood  estimation.  Although  individual  components  were  not 
statistically  significant  and  had  only  little  effects  on  likelihoods 
(with  the  exception  of  eIF4E),  a  combination  gave  improved 
prognostic  power.  This  variable  termed  V  can  be  described  as 
X-Bl/4  +  PB1/2-B2/4,  where  B1  and  PB1  represent  4E-BP1  and 
p4E-BPl.  Each  increase  of  1  in  z  gave  HRs  of  1.15  (P  =  0.006),  1.26 
(P  =  0.002)  and  1.28  (P  =  0.0008)  for  OS,  DFS  and  DSS,  respec¬ 
tively.  This  variable  has  a  highly  significant  relationship  with 
survival  but  this  should  be  treated  with  caution  because  the 
constants  were  determined  using  regressions  for  OS,  and,  thus  to 
an  extent,  significance  is  self  fulfilling,  at  least  for  OS.  The  utility 
of  z,  however,  is  supported  by  the  fact  that  its  relationship  with 
survival  is  more  significant  with  DFS  and  DSS  than  OS,  a  result 
not  predetermined  by  the  approach.  We  also  examined  whether 
y  or  z  give  prognostic  insights  independently  of  NPI  using 
multivariate  analyses.  NPI  and  either  y  or  z  remain  significant 
in  models  for  DFS  (NPI  P <  0.00001;  y  P  =  0.04  or  z  P  =  0.03)  and 
DSS  (NPI  P<  0.00001;  y  P  =  0.02  or  z  P  =  0.02). 

The  statistical  significance  of  relationships  of  y  and  z  with 
survival  show  additional  prognostic  value  from  examining  multi¬ 
ple  eIF4E  pathway  components.  In  addition,  we  have  shown  the 
value  of  these  variables  using  Kaplan-Meier  analyses.  First,  we 
focused  on  patients  with  high  eIF4E  scores  (6  or  7),  as  it  is  in  this 
context  that  differential  expression  of  4E-BPs  would  be  most 
relevant.  Patients  with  eIF4E  scores  6  or  7  have  a  relatively  poor 
prognosis  (Figure  2A-C),  but  no  difference  was  detected  between 
groups  scored  as  6  or  7  in  terms  of  DFS  (Figures  2B  and  3A).  When 
y  was  applied  to  this  cohort  some  discrimination  occurred  with 
improved  prognosis  for  patients  whose  y  scores  were  lowered  by 
4E-BP2  (Figure  3B),  although  the  discrimination  remained 
statistically  nonsignificant.  When  z  was  applied  to  this  cohort 
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Figure  3  Additional  prognostic  information  is  gained  by  combining  assessment  of  eukaryotic  translation  initiation  factor  4E  (elF4E)-binding  proteins 
(4E-BPs)  with  elF4E  analysis.  Kaplan-Meier  survival  analyses  of  disease-free  survival  for  patients  with  either  high  (6  or  7)  elF4E  scores  (A-C)  or  high  (>5.4) 
Nottingham  Prognostic  Index  scores  (D-F).  Patients  were  dichotomised  using  either  elF4E  expression  (A  and  D),  the  y  function  (which  includes  elF4E 
and  4E-BP2  scores)  (B  and  E),  or  the  z  function  (which  includes  scores  for  all  4  markers)  (C  and  F). 


further  discrimination  occurred  (Figure  3C)  showing  how  4E-BPs 
affect  patient  outcome  through  eIF4E.  Second,  we  have  focused  on 
patients  with  high  NPI  (and  consequently  poor  DFS,  Supplemen¬ 
tary  Figure  S6).  These  patients  were  further  stratified  according  to 
eIF4E  expression  (cutoff  5.5  as  suggested  by  the  distribution  in 
Figure  2B)  into  separate  groups  (Figure  3D).  As  before,  when  z  was 
applied  (Figure  3F)  further  discrimination  occurred  allowing 
identification  of  patients  with  very  poor  (high  z),  or  relatively 
good  prognosis  (low  z).  In  this  case,  z  discriminated  into 
statistically  significantly  different  groups  (Figure  3F  Log  rank 
P  =  0.039)  when  use  of  eIF4E  alone  was  not  significant  (Figure  3D 
Log  rank  P  =  0.15).  In  this  casey  was  substantially  less  successful 
as  a  prognostic  indicator  (Figure  3E  Log  rank  P  =  0.5)  showing  the 
importance  of  combining  all  four  components. 


Expressions  of  eIF4E  and  4E-BPs  are  positively  associated 

Although  including  assessments  of  4E-BPs  provided  additional 
prognostic  information  over  that  from  only  eIF4E,  we  were 
surprised  that  influences  of  4E-BPs,  especially  of  4E-BP1,  were 
relatively  weak.  One  explanation  for  this  was  that  expression  levels 
of  eIF4E  and  4E-BPs  were  not  independent  of  each  other.  This 
would  mean  that  that  at  a  given  eIF4E  level,  differential  expression 
of  4E-BPs  -  therefore  differential  modification  of  eIF4E  activity 
-  would  be  relatively  rare,  thereby  minimising  apparent  influences 
of  4E-BPs  in  our  analyses.  Associations  between  marker  expres¬ 
sions  were  examined  using  Spearman’s  correlation  tests  (Table  1). 
eIF4E  expression  showed  moderate  positive  associations  with 
expression  of  both  4E-BPs  (P<  0.0001),  although  weak  positive/no 
association  with  p4E-BPl.  Expression  of  4E-BPs  was  also  mode¬ 
rately  positively  associated  with  each  other  (P  =  0.02).  Expression 
of  p4E-BPl  was  positively  associated  with  4E-BP1  (P<  0.0001) 
(expected  as  4E-BP1  must  be  expressed  to  be  phosphorylated). 


Table  1 

associated 

Expressions 

of  elF4E,  4E-BPI 

and  4E-BP2 

are  positively 

elF4E 

4E-BPI 

4E-BP2 

p4E-BP  1 

elF4E 

_ 

0.31 

0.34 

0.21 

4E-BPI 

0.31 

— 

0.37* 

0.36 

4E-BP2 

0.34 

0.37* 

— 

0.14 

P4E-BPI 

0.21 

0.36 

0.14 

— 

Abbreviations:  el F4E  =  eukaryotic  translation  initiation  factor  4E;  4E-BP  =  elF4E- 
binding  protein;  P4E-BPI  =  phosphorylated  4E-BPI.  Spearman’s  p  correlation  figures 
are  shown.  Associations  are  either  moderately  (*P  =  0.02)  or  highly  significant  (others 
P<  0.0001). 


eIF4E  activity  scores  are  potential  biomarkers  for 
eIF4E-targeted  therapies 

The  eIF4E  pathway  is  a  target  for  cancer  therapy  with  drugs 
that  either  inhibit  eIF4E  activity  indirectly  (mTOR  inhibitors 
that  reduce  4E-BP1  phosphorylation  thereby  inhibiting  eIF4E 
(Smolewski,  2006)),  or  directly  (by  binding  to/reducing  expression 
of  eIF4E  (Graff  et  al>  2008)).  However,  a  concern  with  these  agents 
is  toxicity  resulting  from  general  translational  repression.  Conse¬ 
quently,  selection  of  patients  who  are  most  likely  to  benefit  from 
such  agents  using  predictive  biomarkers  may  aid  their  efficacy. 
Three  selection  criteria  are  apparent:  individuals  with  highest 
eIF4E  levels,  highest  eIF4E  activities  or  highest  levels  of 
phosphorylated  4E-BP1  (especially  relevant  for  mTOR  inhibitors 
as  these  act  by  reducing  4E-BP  phosphorylation).  We  have 
compared  selection  of  potential  treatment  groups  using  these 
criteria  from  our  cohort  for  whom  scores  of  all  four  antigens  were 
available:  first  (group  1),  those  with  eIF4E  scores  of  7  (37/282; 
13.1%);  second  (group  2),  those  with  high  estimated  pathway 
activities  (z^5.75,  cutoff  chosen  to  give  a  similar-sized  group, 
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Figure  4  Use  of  biomarkers  for  dividing  patients  into  potential  treatment 
groups  for  eukaryotic  translation  initiation  factor  4E  (e  I F4E) -directed 
therapy;  different  markers  select  very  different  groups.  A  Venn  diagram 
demonstrating  relationships  between  potential  treatment  groups  selected 
on  the  basis  of  high  elF4E  expression  (I),  high  estimated  elF4E  activity  (‘z’) 
(2),  and  high  p4E-BPI  expression  (3). 

38/282;  13.5%);  finally  (group  3),  those  with  p4E-BPl  scores  of  6  or 
7  (26/282;  9.2%).  We  have  examined  whether  biomarker  choice 
leads  to  selection  of  different  potential  treatment  groups  (Figure  4). 
Substantially  different  patient  groups  were  selected  using  each 
potential  predictive  biomarker.  Only  10  patients  appeared  in  all 
groups  (16%  of  total  patients  selected  in  any  group).  Group  1 
contained  the  largest  proportion  of  uniquely  selected  individuals 
(56%  of  total),  whereas  the  majorities  of  groups  2  and  3  (74  and 
81%  respectively)  overlapped  with  at  least  one  other  group.  Group 
2  (high  z  values)  had  large  overlaps  with  both  other  groups  (46%  of 
group  1  and  81%  of  group  3).  In  addition,  group  2  included  all 
individuals  with  high  eIF4E  and  p4E-BPl  (i.e.  individuals  likely  to 
have  high  eIF4E  activity  by  all  measures)  reflecting  the  fact  that  z 
successfully  takes  account  of  both  eIF4E  and  p4E-BPl,  thereby 
supporting  its  utility  as  a  potential  predictive  marker. 


DISCUSSION 

Expression  of  eIF4E  in  cancer  has  been  studied  extensively, 
however,  expression  does  not  equate  to  activity;  therefore, 
interpretation  of  its  influence  is  more  complex  than  simply 
assessing  expression.  Our  hypothesis  was  that  combined  examina¬ 
tion  of  eIF4E  and  its  regulators  would  allow  greater  insights  into 
eIF4E’s  influence  on  cancer.  Therefore,  we  determined  the 
expression  levels  of  eIF4E  and  its  most  well-established  regulatory 
proteins  4E-BP1,  4E-BP2  and  p4E-BPl  within  tumour  cells  of  a 
large  cohort  of  cancer  patients,  and  have  combined  these  data  into 
an  improved  measure  of  prognosis  and  estimate  of  eIF4E  activity. 
In  common  with  initial  publications  on  eIF4E’s  role  in  cancer 
(Kerekatte  et  al ,  1995;  Li  et  al ,  1998)  and  much  of  the  subsequent 
literature,  we  have  focused  on  breast  cancer. 

There  is  a  wealth  of  literature  on  eIF4E  expression  in  cancer  and, 
although  patient  cohorts  used  have  been  relatively  small  (<200 
individuals),  the  conclusion  that  high  eIF4E  levels  are  associated 
with  poor  prognosis  is  well  established  (De  Benedetti  and  Graff, 
2004).  We  have  also  found  this  (Figure  2A-C)  using  the  largest 
cohort  to  date.  The  4E-BPs  have  received  less  attention.  4E-BP1  is 
more  highly  expressed  in  tumours  than  normal  tissues  (Salehi  and 
Mashayekhi,  2006),  and  expression  correlates  inversely  with 
tumour  progression  (Martin  et  al ,  2000)  -  observations  in 
agreement  with  our  findings.  Despite  the  inverse  correlation  with 
progression,  we  found  examination  of  4E-BP1  expression  not  to 
give  significant  prognostic  insights  (Figure  2D,  Supplementary 
Figure  S5A).  To  our  knowledge,  there  are  no  published  studies 
concerning  4E-BP2  in  cancer  samples.  We  found  that  4E-BP2 
showed  a  nonsignificant  trend  for  association  with  good  prognosis 
(Figure  2E,  Supplementary  Figure  S5B).  High  p4E-BPl  levels  have 


earlier  been  shown  to  correlate  with  grade  and  poor  prognosis 
(Zhou  et  al,  2004;  Castellvi  et  al ,  2006;  Rojo  et  al ,  2007).  In 
contrast,  we  found  detectable  p4E-BPl  (Thr37/46)  in  only  37.7%  of 
the  patients  compared  with  >59%  in  other  reports  (Zhou  et  al , 
2004;  Rojo  et  aly  2007;  Akcakanat  et  aly  2008)),  and  to  provide  little 
prognostic  power  when  analysed  alone  (Figure  2F,  Supplementary 
Figure  S5C).  We  have  carried  out  additional  analyses  with  the  same 
antibody  as  reported  earlier  for  p4E-BPl  (Ser65)  (Zhou  et  aly  2004) 
and  similarly  found  low  expression  (50%  undetectable,  28%  the 
lowest  positive  score).  These  differences  may  relate  to  the  cohorts 
used  as  earlier  studies  had  higher  proportions  of  node  positive 
(Zhou  et  aly  2004),  subsequently  metastatic  (Akcakanat  et  aly  2008) 
or  high-grade  cases  (Rojo  et  aly  2007). 

There  are  very  few  studies  where  multiple  eIF4E  pathway 
components  have  been  analysed.  The  ratios  of  p4E-BPl  to  total  4E- 
BP1  and  of  eIF4E  to  4E-BP1  have  been  shown  to  correlate  with 
high  tumour  grade  (Salehi  and  Mashayekhi,  2006;  Armengol  et  aly 

2007) ,  but  their  relationships  with  survival  were  not  examined.  A 
positive  correlation  between  expression  of  eIF4E  and  p4E-BPl  has 
also  been  noted,  with  the  conclusion  being  that  eIF4E  was  ‘active’ 
in  these  cells  (Nathan  et  aly  2004).  We  have  determined  expression 
levels  of  eIF4E,  4E-BP1,  4E-BP2,  p4E-BPl  and  have  undertaken 
analyses  to  relate  these  together  and  to  survival.  First,  we  found  no 
statistically  significant  relationship  to  survival  for  any  combination 
without  including  eIF4E  -  suggesting  that  eIF4E  is  the  critical 
effector.  Second,  in  support  of  our  initial  hypothesis,  we  showed 
that  combined  assessment  of  the  four  components  allowed 
improved  prognostic  insights  over  eIF4E  alone  (‘z’  function;  text, 
Figure  3).  It  is  important  to  note  that  we  did  not  predetermine 
relationships  between  components  in  z,  rather  these  were  defined 
by  the  best  fit  with  the  data;  the  fact  that  the  relationships  reflect 
our  expectations  from  understanding  the  pathway  supports  the 
view  that  z  is  a  true  estimate  of  eIF4E  activity.  Third,  we  showed 
differential  expression  of  4E-BP2  in  cancer  to  be  more  influential 
in  terms  of  survival  than  4E-BP1.  This  was  shown  by  the 
observations  that  expression  of  4E-BP2,  but  not  4E-BP1,  showed 
a  trend  towards  being  a  prognostic  factor  alone  (Figure  2E, 
Supplementary  Figure  S5B),  provided  an  improved  prognostic 
indicator  in  combination  with  eIF4E  (y),  and  was  the  most 
statistically  significant  component  of  z  after  eIF4E  itself.  This 
observation  may  relate  to  the  fact  that  4E-BP2  binds,  and  therefore 
inhibits  eIF4E  more  strongly  than  4E-BP1  (Abiko  et  aly  2007). 
Interestingly,  we  found  that  expressions  of  eIF4E  and  4E-BPs  were 
positively  associated  (Table  1):  an  unexpected  finding  as  they  are 
functionally  opposed  and  correlate  oppositely  with  grade.  One 
explanation  is  that  4E-BP  translation  may  be  specifically  dere- 
pressed  by  eIF4E’s  action  on  the  5'UTRs  of  their  transcripts, 
representing  a  negative  feedback  loop  within  the  eIF4E  pathway. 

Clinical  trials  of  the  efficacy  and  safety  of  cancer  therapeutics  that 
target  eIF4E  have  been  carried  out  (Graff  et  aly  2008)  and  some 
toxicity  has  been  reported  (O’Donnell  et  aly  2008;  Tabernero  et  aly 

2008) .  Selection  of  individuals  who  are  most  likely  to  benefit  from  the 
agents  may  be  appropriate  in  order  to  avoid  potentially  harmful  and / 
or  ineffective  therapy  in  some  patients.  We  show  that  substantially 
different  patient  groups  are  chosen  using  three  potential  predictive 
biomarkers,  and  therefore  that  use  of  the  best  biomarker  is  important 
for  targeting  of  these  therapies  (Figure  4).  Patients  with  high 
estimated  eIF4E  activity  (‘z’)  (group  2)  and  high  p4E-BPl  levels 
(group  3)  should  provide  good  candidates  for  treatment.  The  former 
group  has  a  particularly  poor  prognosis  (Figure  3C  and  F),  and 
therefore  great  potential  for  clinical  benefit  from  these  drugs. 
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Expression  of  oestrogen  receptor  p  isoforms  is  regulated  by  transcriptional 
and  post-transcriptional  mechanisms 
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Although  ERs  (oestrogen  receptors)  mediate  breast  tumour 
behaviour,  the  precise  role  of  ER/3  remains  unclear.  This  is 
mainly  because  analyses  have  been  complicated  by  the  presence 
in  breast  tissue  of  three  ER/3  protein  variants  (ER/31,  ER/32 
and  ER/35)  that  derive  from  differential  3'  splicing.  We  have 
recently  identified  the  first  known  mechanisms  responsible  for 
the  differential  control  of  isoform  expression,  involving  regulation 
of  translation  via  5 -UTRs  (untranslated  regions).  In  the  present 
study,  we  have  uncovered  further  complexity  involving  the 
influence  of  multiple  promoters  and  cross-talk  between  5'-  and 
3  -UTRs.  We  demonstrate  that  full-length  ER/3  mRNAs  are  trans¬ 
cribed  from  three  separate  promoters;  two  promoters  are  well- 
established  within  the  literature,  whereas  the  third  represents  a 
novel  finding.  Each  promoter  produces  transcripts  with  distinct  5'- 
UTRs.  The  differential  3'  splicing  that  produces  transcripts  coding 


for  the  ER/3  isoforms  also  defines  isoform-specific  3-UTRs.  We 
identified  exact  3/-UTR  sequences  for  each  isoform,  and  have 
shown  that  alternative  polyadenylation  sites  are  used  in  a  cell-type 
specific  manner  to  produce  transcripts  with  3/-UTRs  of  different 
lengths.  Critically,  we  show  that  5'-  and  3/-UTRs  combine  to 
specify  the  efficiencies  with  which  individual  transcripts  are 
translated,  with  3/-UTR  length  having  a  key  influence.  In  addition, 
we  demonstrate  how  17/3-oestradiol,  a  key  driver  of  breast  cancer 
development,  affects  the  regulation  of  ER/3  expression  at  both 
transcriptional  and  translational  levels. 


Key  words:  breast  tissue,  differential  polyadenylation,  esr2 , 
oestrogen  receptor  (ER),  translational  regulation,  untranslated 
region  (UTR). 


INTRODUCTION 

ERs  (oestrogen  receptors)  are  critical  mediators  of  oestrogen 
function  and  play  roles  in  many  pathological  processes,  especially 
carcinogenesis  [1].  Although  the  roles  of  ERa  are  relatively  well 
understood,  those  of  ER/3  remain  unclear.  This  is  due,  in  part,  to 
a  reported  discrepancy  between  ER/3  expression  at  mRNA  and 
protein  levels,  leading  to  conflicting  expression  data  [2,3].  Also, 
ER/3  is  expressed  as  at  least  five  protein  isoforms,  derived  from 
differential  3'  splicing  of  ER/3  transcripts  [4],  yet  their  potential 
to  have  distinct  functions  has  often  been  ignored  in  favour  of 
analysis  of  total  ER/3.  In  breast  tissue,  ER/31,  ER/32  and  ER/35 
predominate,  with  each  appearing  to  have  separate  biological 
functions  as  demonstrated  by  their  associations  with  different 
breast  cancer  types  and  prognoses  [5-8].  Differences  between  the 
isoforms  are  also  evident  in  terms  of  comparisons  of  expression 
in  normal  breast  tissue  and  breast  cancer.  ER/31  is  frequently 
down-regulated  in  cancer  compared  with  normal  cells  [9,10], 
suggesting  that  it  may  function  as  a  tumour  suppressor  [11-13]. 
However,  ER/32  appears  to  be  up-regulated  during  carcinogenesis 
[10,14],  whereas  ER/35  may  also  be  up-regulated,  at  least  at 
the  level  of  mRNA  [15].  Little  is  known  about  the  mechanisms 
responsible  for  these  changes  in  ER/3  expression.  ER/3  promoter 
methylation  has  frequently  been  observed  in  breast  cancers,  and 
this  is  thought  to  be  responsible  for  down-regulation  of  some  ER/3 
transcripts,  although  this  appears  to  be  at  odds  with  the  reported 
up-regulation  of  ER/32  and  ER/35.  We  have  recently  identified 
the  first  known  mechanisms  responsible  for  differential  control 
of  expression  of  the  different  isoforms  involving  regulation  of 
translation  via  two  alternative  5'-UTRs  (untranslated  regions) 
[16].  We  have  now  extended  this  theme  to  examine  thoroughly 


the  regulatory  functions  of  the  extensive  range  of  ER/3  UTRs,  and 
have  determined  further  mechanisms  for  differential  control  of 
ER/3  isoform  expression. 

EXPERIMENTAL 

Cell  culture,  transfection,  flow  cytometry  and  dual  luciferase 
assays 

Cell  lines  representing  breast  cancers  of  luminal  (MCF7)  and 
basal  (MDAMB-231)  subtypes  and  benign  non-transformed 
breast  tissue  (HB2)  were  obtained  from  the  European  Collection 
of  Animal  Cell  Cultures.  Cells  were  cultured  in  RPMI  1640 
medium  containing  5  %  (v/v)  FBS  (fetal  bovine  serum)  (MCF7, 
MDAMB-23 1)  or  DMEM  (Dulbecco’s  modified  Eagle’s  medium) 
containing  or  10%  (v/v)  FBS  (HB2)  (both  Invitrogen)  at  37  °C 
under  5%  C02.  Bi-monthly  mycoplasma  checks  (Myco Alert® 
Mycoplasma  detection  assay,  Lonza)  were  consistently  negative 
and  STR  (short  tandem  repeat)  profiles  confirmed  cell  identity. 
Cells  were  transfected  as  described  previously  [16].  For 
experiments  using  exogenous  E2  (17^-oestradiol),  cells  were 
cultured  in  Phenol  Red-free  medium  (Invitrogen)  supplemented 
with  charcoal- stripped  FBS.  At  5  h  post- transfection,  fresh 
medium  or  medium  containing  10  nM  E2  (Sigma- Aldrich)  was 
added  before  analysis  after  24  h.  For  flow  cytometry,  cells  were 
removed  from  wells  with  trypsin  and  resuspended  in  fresh  medium 
containing  1  %  (v/v)  serum.  GFP  (green  fluorescent  protein) 
expression  was  quantified  (mean  fluorescent  intensity  of  104 
events  after  exclusion  of  debris/dead  cells  on  the  basis  of  forward 
activated  light  scatter  against  side  scatter)  at  525  nm  (LSRII, 
BD  Biosciences).  Gates  were  set  so  that  <1  %  of  untransfected 


Abbreviations  used:  E2,  17/J-oestradiol;  ER,  oestrogen  receptor;  FBS,  fetal  bovine  serum;  GFP,  green  fluorescent  protein;  ORF,  open  reading  frame; 
qPCR,  quantitative  PCR;  RACE,  rapid  amplification  of  cDNA  ends;  RT,  reverse  transcriptase;  uORF,  upstream  ORF;  UTR,  untranslated  region. 

1  To  whom  correspondence  should  be  addressed  (email  t.hughes@leeds.ac.uk). 
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cells  were  defined  as  expressing  GFP.  Dual-luciferase  assays 
(Promeg a)  were  performed  according  to  the  manufacturer’s 
instructions  using  a  Lumat  LB 95 07  luminometer  (Berthold 
Technologies).  pSV40 -Renilla  (Promega)  was  used  as  a  control. 

Plasmid  construction 

pTH-GFPa  and  GFP  reporters  for  UTRa  and  UTRc  have  been 
described  previously  [16,17].  The  GFP  reporter  for  the  El 
5-UTR  was  cloned  similarly  to  those  for  UTRa  and  UTRc; 
the  UTR  was  amplified  by  PCR  from  cDNA  prepared  from 
MCF7  cells  and  was  cloned  upstream  of  the  GFP  ORF  (open 
reading  frame)  in  pTH-GFPa.  Importantly,  this  strategy  removes 
the  5'  end  of  the  multiple  cloning  site,  allowing  the  inserted 
5  -UTRs  to  be  immediately  adjacent  to  the  transcriptional  start 
site.  3-UTRs  were  amplified  by  PCR  from  cDNA  prepared  from 
HB2  cells  and  were  cloned  downstream  of  the  GFP  ORF  in 
each  of  the  5 '-UTR  reporter  constructs  with  BamHI/Hindlll. 
Promoter  sequences  were  amplified  from  MCF7  genomic  DNA 
and  cloned  into  pGL3 -Basic  (Promega)  using  KpnI/Nhel.  The 
El  promoter  was  cloned  as  two  fragments:  a  3'  fragment 
cloned  using  KpnI/Nhel,  then  a  5'  fragment  cloned  using  KpnI. 
Primer  sequences  are  listed  in  Supplementary  Table  SI  at 
http://www.BiochemJ.org/bj/429/bj4290283add.htm. 

cDNA  synthesis  and  PCR 

RNA  was  purified  from  cells  with  RNeasy  kits  (Qiagen); 
contaminating  DNA  was  removed  with  Turbo  DNase  I 
(Applied  Biosystems).  First-strand  cDNA  was  synthesized 
using  Superscript  II  (according  to  the  manufacturer’s  protocol) 
and  oligo(dT)  or  random  hexamers.  Triplicate  real-time  PCR 
analysis  was  performed  (Applied  Biosystems  SYBR®  Green 
PCR  Master  Mix  and  7900HT  machine).  Dissociation  curves 
and  serial  cDNA  dilutions  were  performed  to  ensure  primer 
specificities  and  equivalent  amplification  efficiencies;  correlation 
coefficients  of  >0.985  and  primer  efficiencies  of  >95%  and 
<  100  %  were  deemed  acceptable.  Reactions  were  also  performed 
using  template  lacking  RT  (reverse  transcriptase):  products 
were  either  undetectable  or  greatly  reduced  (> 30000-fold  less 
product  than  the  equivalent  RT+),  hence  genomic  or  plasmid 
DNA  contamination  was  not  considered  to  interfere  with  data. 
Expression  of  UTRs  was  determined  relative  to  expression  of  the 
RPLP0  (36B4)  gene  [18].  RACE  (rapid  amplification  of  cDNA 
ends)  was  performed  using  5  RACE  System  2  or  3'RACE  System 
(both  Invitrogen)  according  to  the  manufacturer’s  instructions 
(both  standard  and  a  modified  adapter  primer  were  used  for  ER/3  2 
3'RACE).  All  primer  sequences  are  listed  in  Supplementary  Table 
SI.  Products  were  analysed  on  2.5%  agarose  gels  [0.5  /xg/ml 
ethidium  bromide,  1  x  TBE  (45  mM  Tris/borate  and  1  mM 
EDTA)]  and  visualized  on  an  UV  transilluminator.  Products  were 
excised  from  gels  and  cloned  into  pGEM-Teasy  (Promega);  at 
least  five  clones  for  each  were  sequenced.  Note  that  products 
are  larger  than  the  UTRs  they  represent  since  they  include  some 
reading  frame  and  the  RACE  adapters. 

RNA  structure  and  statistical  analyses 

Modelling  was  performed  using  mfold  v3.1  to  predict  potential 
secondary  structures  for  RNA  molecules  as  described  previously 
[16].  The  algorithm  finds  base-pairing  solutions  that  are  sterically 
possible  and  release  the  greatest  amount  of  free  energy  (AG) 
during  structural  folding;  more  stable  structures  release  more 
energy  as  they  form  and  therefore  have  greater  AG  values  [19]. 
Student’s  t  test  was  used  for  statistical  analysis  using  Microsoft 


Excel.  All  P  values  were  two-sided;  P  <  0.05  was  considered 
significant. 

RESULTS 

A  novel  5'-UTR  for  ER/3 

We  studied  previously  the  regulatory  roles  of  two  ER/3  5 -UTRs 
that  we  termed  UTRa  and  UTRc  [16].  These  5-UTRs  result  from 
transcription  initiation  from  two  alternative  promoters  and  mutu¬ 
ally  exclusive  splicing  of  the  untranslated  first  exons,  exon  OK  or 
ON,  to  the  first  coding  exon  (exon  1)  [20].  5 -RACE  analyses  were 
performed  in  breast  cell  lines  in  order  to  examine  ER/3  5'-UTRs 
further.  Products  representing  5  -UTRs  containing  only  sequence 
from  exon  1  or  from  immediately  upstream  of  the  5'  splice  site 
of  exon  1  were  the  sole  products  obtained  from  primers  located 
within  the  coding  region  of  exon  1  (Figure  1  A).  UTRa  or  UTRc 
could  only  be  detected  by  5-RACE  using  primers  specific  for 
their  respective  upstream  exons,  OK  and  ON  [16].  RACE  products 
containing  only  sequence  from  within  the  accepted  bounds  of  exon 
1  could  be  interpreted  in  two  ways:  as  truncated  versions  of  UTRa 
or  UTRc,  or  as  complete  5'-UTRs  derived  from  transcriptional 
initiation  within  exon  1 .  Products  containing  sequence  upstream 
of  the  accepted  bounds  of  exon  1  were  likely  to  result  from 
transcriptional  initiation  adjacent  to  exon  1.  Figure  1(B)  shows 
an  alignment  of  mRNAs  containing  UTRa,  UTRc  or  the  novel 
shorter  5-UTR  (‘UTR-E1’)  with  the  5'  end  of  the  human  ER/3 
gene  ( esr2 )  on  14q23.  The  diagram  depicts  putative  transcriptional 
initiation  (bent  black  arrows)  over  a  range  of  sequences  upstream 
of  and  within  exon  1  allowing  the  expression  of  a  5'  extension 
to  exon  1  (grey  box)  on  some  UTR-E1 -containing  transcripts. 
We  next  performed  qPCR  (quantitative  PCR)  to  examine  the 
relative  expression  of  UTRa,  UTRc  and  UTR-E1  in  breast  cell 
lines.  We  analysed  expression  of  UTR-E1  using  a  primer  that 
was  complementary  to  the  sequence  within  this  5'  extension  of 
exon  1,  and  was  therefore  not  contained  within  UTRa  or  UTRc. 
It  is  worth  noting  that  these  analyses  may  under-represent  UTR- 
E1  expression,  since  RACE  analyses  show  that  the  5'  extension 
of  exon  1  is  present  on  only  a  subset  of  UTR-E1 -containing 
transcripts.  All  three  breast  cell  lines  examined  expressed  UTR-E1 
(Figure  1C).  Products  were  not  amplified  from  mock  reverse  tran¬ 
scription  reactions,  providing  validation  that  transcribed/reverse- 
transcribed  sequences  were  detected  rather  than  contaminating 
genomic  DNA.  In  HB2  and  MCF7  cells,  UTRa  was  the  majority 
species,  with  UTRc  and  UTR-E1  being  expressed  at  similar 
lesser  levels.  In  contrast,  MDAMB-231  cells  expressed  ~  12.5- 
fold  more  UTRc,  and  ~  3.5-fold  more  UTR-E1,  than  UTRa.  We 
concluded  that  UTR-E1  represented  a  novel  5-UTR  for  ER/3. 

A  third  transcriptional  promoter  for  ER/3 

Next,  our  hypothesis  was  that  the  DNA  immediately  upstream 
of  exon  1  acts  as  a  promoter  allowing  expression  of  transcripts 
containing  UTR-E1.  We  cloned  ~2kb  of  the  genomic  DNA 
directly  upstream  of  exons  OK,  ON  or  El  into  separate  luciferase 
reporter  vectors  and  performed  luciferase  assays  in  HB2  and 
MCF7  breast  cell  lines  in  order  to  examine  the  relative  promoter 
activities  of  the  two  known  ER/3  promoters  and  the  putative 
third  promoter.  Activities  are  shown  in  each  cell  type  (Figure  2). 
The  putative  third  promoter  (‘promoter  El’)  showed  activity  in 
both  cell  lines.  In  HB2  cells,  promoter  El  showed  intermediate 
activity  between  promoters  ON,  the  most  active,  and  OK,  the  least 
active.  In  MCF7  cells,  promoters  OK  and  ON  determined  similar 
transcriptional  activities,  whereas  promoter  El  was  ~70% 
weaker.  We  concluded  that  ER/3  transcription  is  indeed  driven 
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Figure  1  A  novel  ER/?  5'-UTR 

(A)  5'-RACE  analyses  of  ER£  transcripts  were  performed;  products  representing  a  previously  unidentified  5'-UTR  containing  sequence  from  or  immediately  upstream  of  exon  1  were  amplified. 
Products  from  two  breast  cell  lines  are  shown  along  with  a  diagram  representing  ER£  exon  1  (genomic  locations  with  respect  to  the  translational  start  are  indicated).  (B)  Alignment  of  the  5'  end  of 
the  human  ER>0  gene  with  mRNAs  containing  either  UTRa,  UTRc  or  the  novel  UTR-E1 .  Coding  regions  (black  boxes),  transcriptional  start  sites  (black  arrows)  and  primers  used  for  the  qPCR  analysis 
of  each  5'-UTR  (black  arrows)  are  shown.  The  grey  box  indicates  sequence  included  in  a  subset  of  transcripts  with  UTR-E1.  (C)  Relative  expression  of  UTRa,  UTRc  and  UTR-E1  in  breast  cell  lines 
was  examined  by  qPCR;  all  three  cell  lines  examined  expressed  UTR-E1. 


promoters 


promoters 


Figure  2  A  third  ER /?  promoter  shows  activity  in  breast  cell  lines 

Approx.  2  kb  of  genomic  DNA  directly  upstream  of  exons  OK,  ON  or  El  was  cloned  into 
luciferase  reporter  vectors.  HB2  and  MCF7  cell  lines  were  transiently  transfected  with  equal 
copy  numbers  of  luciferase  reporter  lacking  additional  promoter  sequences  (empty)  or  containing 
promoter  sequences  as  shown,  and  luciferase  assays  were  performed  after  24  h.  A  minimum 
of  two  independent  experiments  were  performed,  and  within  each  experiment  three  technical 
replicates  were  included.  Results  are  means  ±  S.D.  of  technical  triplicates  within  a  representative 
experiment. 


by  a  third,  previously  uncharacterized,  promoter  that  determines 
the  expression  of  UTR-E1. 


5'-UTRs  differentially  regulate  efficiency  of  ER/?  translation 

We  have  shown  previously  that  UTRa  and  UTRc  have  profound 
and  differential  influences  on  ER/3  translation.  We  hypothesized 
that  UTR-E1  may  also  influence  translation,  thus  we  extended  our 
previous  analyses  to  test  this  using  our  established  GFP  reporter 
assay  [16,21].  UTRa,  UTRc  and  UTR-E1  were  cloned  upstream 
of  the  GFP  reading  frame  in  expression  vectors.  For  UTR-E1,  we 
cloned  the  sequence  encoded  by  the  published  extent  of  exon  1, 
representing  a  commonly  identified  5 -RACE  product.  Cells  were 
transiently  transfected  with  equal  numbers  of  copies  of  vectors 
to  allow  expression  of  GFP  mRNAs  either  with  non-regulatory 
5 -UTRs  (positive  control;  ‘con’),  or  with  UTRa,  UTRc  or  UTR- 
El.  GFP  protein  expression  was  measured  by  flow  cytometry, 
and  GFP  mRNA  expression  was  measured  by  qPCR  allowing 
determination  of  relative  translational  efficiencies  for  each  GFP 
mRNA  (Figure  3).  As  published  previously  [16],  UTRa  and  UTRc 
inhibited  translation,  with  UTRa  being  strikingly  inhibitory  in 
nature.  UTR-E1  was  also  inhibitory  in  nature  despite  being  short 
(only  90  nucleotides;  as  compared  with  UTRa,  289  nucleotides, 
and  UTRc,  418  nucleotides)  and  having  a  very  low  predicted 
degree  of  secondary  structure  [as  assessed  by  the  theoretical 
change  in  free  energy  of  folding,  AG;  UTR-El’s  predicted  AG  is 
— 14  kcal/mol  (1  kcal  =  4. 184  kJ),  whereas,  for  comparison,  those 
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5’UTR:  con  UTRa  UTRc  UTR-E1 


Figure  3  5'-UTRs  regulate  ER>S  translational  efficiency 

Reporters  were  constructed  to  express  mRNAs  containing  the  GFP  reading  frame  preceded  by 
different  S'-UTRs:  a  control  (con)  sequence  lacking  regulatory  motifs,  UTRa,  UTRc  or  UTR-E1. 
Cell  lines  were  transiently  transfected  with  equal  copy  numbers  of  either  control  or  experimental 
constructs.  GFP  protein  and  mRNA  were  quantified  by  flow  cytometry  and  real-time  PCR 
respectively.  Translational  efficiency  (protein  synthesized  per  unit  of  mRNA)  is  presented  relative 
to  the  GFP  control  that  lacks  a  specialized  5'-UTR.  A  minimum  of  two  independent  experiments 
were  performed  and  within  each  experiment  three  technical  replicates  were  included.  Results 
are  means  ±  S.D.  of  technical  triplicates  within  a  representative  experiment. 

of  UTRa  and  UTRc  are  —84  and  —166  kcal/mol].  A  consistent 
pattern  of  relative  influences  for  each  UTR  was  seen  in  all  cell 
lines,  with  UTRa  being  more,  and  UTR-E1  less,  inhibitory.  We 
concluded  that  ER/3  5-UTRs  specified  the  efficiencies  with  which 
downstream  ORFs  are  translated.  We  also  noted  that  promoter  El 
may  be  especially  important  in  terms  of  defining  ER/3  function, 
since  transcripts  from  this  promoter,  which  contain  UTR-E1,  are 
the  most  efficiently  translated  ER/3  mRNAs. 


Cross-talk  between  5'-  and  3'-UTRs  influences  translation  of 
individual  ER/3  isoforms 

Transcripts  for  at  least  three  functionally  distinct  ER/3  isoforms 
are  produced  in  breast  cells  [15].  These  are  derived  from 
differential  3'  splicing  of  the  final  ER/3  exons  [4].  It  is  evident 
that  this  differential  splicing  must  also  confer  different  3  -UTRs 
on  transcripts  for  each  isoform,  although  the  exact  3 -UTR 
sequences  are  poorly  defined,  especially  with  respect  to  which 
potential  polyadenylation  sites  are  used.  We  were  interested  to 
examine  whether  ER/3  3 -UTRs  might  influence  translation.  First, 
we  performed  3' -RACE  analyses  for  each  isoform  in  HB2  and 
MCF7  breast  cells  to  identify  3-UTR  sequences  (Figure  4A). 
For  ER/31,  we  were  only  able  to  amplify  a  product  in  MCF7 
cells,  probably  because  ER/3 1  expression  levels  were  low  in  HB2 
cells.  The  3-UTR  identified  was  of  242  nucleotides,  representing 
a  considerable  3'  extension  of  the  published  sequence  (108 
nucleotides;  GenBank®  accession  NM_001437.2).  For  ER/32, 
we  identified  three  3-UTRs  of  similar  lengths:  an  85  nucleotide 
sequence  in  HB2  cells  and  103  or  108  nucleotide  sequences 
in  MCF7  cells;  each  apparently  terminating  at  polyadenylation 
sites  slightly  more  proximal  than  the  published  3-UTR  sequence 
of  120  nucleotides  (GenBank®  accession  number  AF051428). 
Cell- type- specific  alternative  3-UTRs  were  identified  for  ER/3 5: 
sequences  of  234  nucleotides  in  HB2  and  of  79  nucleotides  in 
MCF7  cells.  These  represented  the  use  of  either  proximal  or  distal 
polyadenylation  sites  as  compared  with  the  published  sequence  of 
177  nucleotides  (GenBank®  accession  numbers  DQ838583.1  and 
AF061055.1).  We  did  not  detect  expression  of  an  extended  ER/32 
3-UTR  that  is  represented  within  GenBank®  (accession  number 
NM_001040276.1;  498  nucleotides).  However,  RACE  reactions 
can  be  biased  towards  amplification  of  shorter  sequences, 
therefore  we  designed  primers  to  assess  expression  of  this  specific 
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Figure  4  Differential  splicing  confers  multiple  and  different  3'-UTRs  on  transcripts  for  each  ER (3  isoform 

(A)  3'-RACE  analyses  for  each  ER£  isoform  were  performed  on  the  breast  cell  lines  HB2  (FI)  and  MCF7  (M).  Products  were  identified  and  sequenced.  PCR  analysis  was  also  performed  using  primers 
specific  for  an  ER^2  3'-UTR  that  is  represented  within  GenBank®  (/32-long;  right-hand  panel).  (B)  Alignment  of  the  3'  end  of  the  human  ERyB  gene  with  mRNAs  for  each  isoform  containing  the 
3'-UTRs  (open  boxes)  identified.  Coding  regions  (black  boxes)  and  the  location  of  polyadenylation  site  [p(A)]  and  primers  used  for  amplification  of  products  in  (A)  (black  arrows)  are  shown. 
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Figure  5  Isoform-specific  3'-UTRs  define  different  translational  efficiencies  in  conjunction  with  5'-UTRs 

S'-UTRs  were  cloned  downstream  of  the  GFP  reading  frame  in  each  5'-UTR  GFP  reporter  construct.  HB2  (A)  and  MCF7  (B)  cells  were  transiently  transfected  with  reporter  constructs  and  translational 
efficiencies  determined  relative  to  the  GFP  control.  A  minimum  of  two  independent  experiments  were  performed  and  within  each  experiment  three  technical  replicates  were  included.  Results  are 
means  ±  S.D.  of  technical  triplicates  within  a  representative  experiment. 


3 -UTR.  The  longer  ER02  3-UTR  (‘02-long’)  was,  in  fact, 
expressed  in  both  cell  lines  (Figure  4A;  right-hand  panel).  These 
data  are  summarized  in  Figure  4(B)  in  the  form  of  an  alignment 
of  the  y  ends  of  mRNAs  for  each  isoform  with  the  3'  end  of  the 
human  ER0  gene. 

We  were  interested  to  examine  whether  these  isoform- specific 
3-UTRs  define  different  translational  efficiencies  in  conjunction 
with  the  5'-UTRs,  thereby  allowing  differential  expression  of  the 
isoforms.  We  cloned  3-UTRs  downstream  of  the  GFP  reading 
frame  in  each  5 -UTR  GFP  reporter  construct.  For  ER02,  we 
focused  on  comparison  of  /32-long,  as  detected  by  specific 
PCR,  with  the  sequences  identified  using  RACE  (‘/32-short’; 
we  examined  a  120  nucleotide  sequence,  GenBank®  accession 
number  AF051428).  For  ER/35,  we  examined  both  sequences 
identified  by  RACE  (‘05  HB’  and  ‘05  MCF’).  HB2  and  MCF7 
cells  were  transiently  transfected  with  equal  numbers  of  copies  of 
vectors  to  allow  expression  of  GFP  mRNAs  with  non-regulatory 
UTRs  (positive  control;  ‘con’),  or  with  reporters  for  each  isoform- 
specific  3-UTR  in  combination  with  each  5-UTR.  Translational 
efficiencies  were  determined  as  described  above  and  are  presented 
relative  to  positive  controls  (Figure  5).  3-UTRs  had  potent  and 
differential  influences  on  translational  efficiencies.  The  ER/31  3'- 
UTR  had  little  influence  on  the  translational  efficiencies  specified 
by  the  5 -UTRs  in  either  cell  line  (compare  lane  2  with  lane  1, 
lane  8  with  lane  7  and  lane  14  with  lane  13).  In  contrast,  the  two 
alternative  ER/32  3-UTRs  had  markedly  different  effects.  02- 
long  had  little  influence  on  the  translational  efficiencies  specified 
by  5-UTRs  (compare  lane  4  with  lane  1,  lane  10  with  lane  7  and 
lane  16  with  lane  13).  However,  02-short  induced  dramatic  and, 
in  the  case  of  MCF7  cells,  total  derepression  of  the  translation 


inhibition  specified  by  5-UTRa  (compare  lane  3  with  lane  1; 
P  <  0.0001  in  both  cell  lines),  but  had  relatively  little  influence 
when  paired  with  UTRc  or  UTR-E1  (compare  lane  9  with  lane  7 
and  lane  15  with  lane  13).  ER05  3-UTRs  had  smaller  influences 
on  translational  efficiencies.  05  MCF  induced  derepression  of  the 
inhibitory  influence  of  both  UTRc  and  UTR-E1  (compare  lane  1 1 
with  lane  7  and  lane  13  with  lane  8;  P  <  0.02  in  all  cases),  whereas 
05  HB  had  little  influence  except  for  inducing  strong  translation 
when  paired  with  UTR-E1  (compare  lane  18  with  lane  13)  in  only 
HB2  cells  (notably,  these  are  the  cells  from  which  this  3-UTR 
was  cloned).  We  concluded  that  cross-talk  between  5'-  and  3'- 
UTRs  had  profound  influences  on  the  translational  efficiency  of 
transcripts  for  individual  ER0  isoforms,  with  length  of  3-UTRs, 
as  defined  by  differential  use  of  polyadenylation  sites,  being  a 
critical  factor  in  determining  the  outcome. 

E2  modifies  ER0  transcription  and  translation 

E2  influences  the  downstream  effects  of  ERs  by  binding  to  and 
modifying  their  activity  as  transcription  factors.  It  has  also  been 
reported  that  E2  has  an  impact  on  expression  of  ER0  itself  at  both 
transcriptional  [22,23]  and  translational  [24]  levels  by  largely 
unknown  mechanisms.  We  next  investigated  the  influence  of 
E2  on  ER0  transcription  and  translation  using  our  reporters  for 
both  of  these  regulatory  stages.  Cells  were  transfected  with  equal 
copy  numbers  of  either  luciferase  reporters  containing  different 
ER0  promoters  (as  in  Figure  2)  or  GFP  reporters  containing 
isoform-specific  UTR  pairings  (as  in  Figure  5)  for  analysis 
of  influences  on  transcription  and  translation  respectively. 
Cells  were  then  treated  with  vehicle  or  E2  for  24  h  and  were 
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B  Promoter  ON 
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C  Promoter  El 


HB2  MCF7 

Figure  6  E2  modifies  ER/3  transcription 

HB2  and  MCF7  cells  were  transfected  with  equal  copy  numbers  of  luciferase  reporters  containing 
different  ER£  promoters  (as  in  Figure  2).  Cells  were  treated  with  vehicle  (white  bars)  or  10  nM  E2 
(black  bars)  for  24  h,  and  luciferase  assays  were  performed.  E2  caused  an  increase  in  the  activity 
of  promoter  OK  in  MCF7  cells,  but  not  in  FIB2  cells  (A),  whereas  it  had  no  significant  effect  on 
promoters  ON  or  El  (B  and  C).  A  minimum  of  two  independent  experiments  were  performed 
and  within  each  experiment  three  technical  replicates  were  included.  Results  are  means  ±  S.D. 
of  technical  triplicates  within  a  representative  experiment. 


analysed  for  luciferase  activity  or  translational  efficiency  of 
GFP  transcripts.  In  terms  of  transcriptional  activity,  E2  caused  a 
1 .6-fold  activation  of  the  OK  promoter  in  MCF7  cells  (P  =  O.OOl), 
but  not  in  HB2  cells  (Figure  6A),  while  having  no  significant 
effect  on  promoters  ON  or  El  (Figures  6B  and  6C).  Effects  of  E2 
on  translational  efficiencies  specified  by  the  various  UTR  pairings 
were  more  complex.  Translation  of  transcripts  containing  the 
ER/31  3-UTR  were  increased  by  E2  when  combined  with 
UTRc  or  UTR-E1  (Figure  7A;  P  =  0.008  and  0.03  respectively). 
ER^2  3 -UTRs  (both  short  and  long)  specified  an  E2-dependent 
decrease  in  translational  efficiency  when  paired  with  5-UTRa 
(Figure  7B;  both  P  =  0.003),  but  an  increase  when  paired  with 
UTR-E1  (P  =  0.01  and  0.049  respectively).  ER^5  3-UTRs  also 
specified  differential  responses  to  E2,  determining  both  decreases 
0 35  MCF  paired  with  UTR-E1;  P  =  0.005)  and  increases  in 


3‘UTR:  con  [31  pi  pi 


~  0 


5fUTR:con  UTRa  UTRc  UTR-E1 

3T  UTR  icon  02s  [321  p2s  02)  02s  021 


3!UTR:con  05M  05H  p5M  p5H  05M  05H 


Figure  7  E2  modifies  ER/?  translation 

MCF7  cells  were  transfected  with  equal  copy  numbers  of  GFP  reporters  containing  UTR  pairings 
(as  in  Figure  5).  Cells  were  treated  with  vehicle  (white  bars)  or  10  nM  E2  (black  bars)  for  24  h, 
and  translational  efficiency  of  GFP  transcripts  was  determined.  A  minimum  of  two  independent 
experiments  were  performed,  and  within  each  experimentthree  technical  replicates  were  included. 
Results  are  means  ±  S.D.  of  technical  triplicates  within  a  representative  experiment. 

translation  (ft 5  HB  paired  with  UTRc;  P  =  0.002).  We 
concluded  that  E2  has  only  relatively  weak  influences  on  ER/3 
transcription  and  translation.  However,  as  these  influences  act 
differentially  on  the  promoters  and  on  the  isoform-specific  UTR 
pairings,  these  weak  influences  cumulatively  have  potential  to 
modify  the  balance  of  ER/3  isoforms. 
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DISCUSSION 

It  is  well  established  that  the  ER/3  gene  has  two  promoters 
allowing  the  expression  of  full-length  transcripts,  promoters  OK 
and  ON  [20].  mRNAs  transcribed  from  each  promoter  have 
different  5 -UTRs,  termed  UTRa  and  UTRc  [16]  (Figure  IB). 
However,  full-length  ER/3  transcripts  containing  neither  of  these 
UTRs  have  been  reported  [25],  and  the  presence  of  an  additional 
promoter,  giving  rise  to  further  alternative  5  -UTRs,  has  long 
been  suspected  [26].  We  have  identified  an  additional  5-UTR, 
UTR-E1,  and  the  corresponding  novel  promoter  region,  promoter 
El  (Figures  1  and  2).  A  number  of  previous  studies  have 
attempted  to  assess  the  relative  transcriptional  strengths  and 
thereby  the  importance  of  promoters  OK  and  ON  [27,28] .  However, 
our  previous  data,  demonstrating  that  ER/3  5-UTRs  determine 
differential  and  cell-type- specific  translational  inhibition  [16], 
cast  doubt  on  these  assessments.  Reporter  assays  including  exonic 
sequences  would  combine  translational  differences  with  the 
expected  measure  of  transcription.  We  have  therefore  investigated 
promoter  activities  of  only  sequences  immediately  upstream  of  the 
exons  thereby  avoiding  translational  regulatory  motifs  (although 
risking  failure  to  include  important  transcription-factor-binding 
sites  within  the  exons),  and  have  analysed  translational  regulation 
separately.  We  found  promoter  ON  to  be  the  most  active  in  both 
cell  lines  tested,  although  promoters  El  and  OK  are  also  active 
(Figure  2).  Methylation  of  CpG  islands  within  promoter  ON  has 
been  reported  in  cancers,  including  those  of  the  breast  [29],  ovary 
[30]  and  prostate  [28],  and  this  is  thought  to  be  responsible  for 
down-regulation  of  some  forms  of  ER/3  during  carcinogenesis 
[9,10,31].  In  this  context,  it  is  interesting  to  note  the  presence 
of  a  CpG  island  within  promoter  El  (from  —2180  to  —1878 
with  respect  to  the  translational  start  site);  our  preliminary  data 
suggest  that  these  sequences  are  also  commonly  methylated  in 
breast  cancers  [32].  It  is  also  interesting  that  the  relative  activities 
of  the  promoters  as  determined  by  luciferase  assay  (Figure  2) 
do  not  correlate  with  the  expression  of  the  5-UTRs  derived 
from  the  endogenous  promoters  (Figure  1C).  This  may  relate 
to  the  inhibitory  methylation  of  endogenous  promoters,  or  to  the 
influence  of  transcription  factors  binding  outside  the  sequence 
included  in  the  reporter  constructs. 

We  demonstrated  previously  that  ER/3  expression  is  regulated 
at  the  level  of  translation,  with  UTRa  and  UTRc  having  potent  and 
differential  influences  [16].  When  these  analyses  were  extended 
to  include  UTR-E1,  we  found  that  UTR-E1  allowed  the  most 
efficient  translation  of  the  ER/3  5-UTRs  (Figure  3).  Therefore 
transcripts  from  promoter  El,  which  contain  UTR-E1,  are  likely 
to  contribute  disproportionally  to  ER/3  protein  levels.  Moreover, 
a  key  finding  of  the  present  study  is  that  translational  regulation 
acts  differentially  on  transcripts  for  the  ER/3  isoforms  1,  2  and  5. 
These  transcripts  are  derived  from  differential  3'  splicing  [4], 
giving  them  isoform- specific  3-UTRs.  We  identified  3 -UTR 
sequences  for  each  isoform  (Figure  4)  and  demonstrated  that, 
when  combined  with  ER/3  5-UTRs,  these  specified  a  wide  range 
of  translational  efficiencies  (Figure  5).  Most  striking  was  the 
ability  of  different  ER/32  3/-UTRs  to  define  up  to  a  10-fold 
difference  in  translational  efficiency  when  paired  with  UTRa.  5'- 
UTRa  inhibits  translation  on  account  of  inefficient  translational 
scanning  induced  by  secondary  structure  and  initiation  at  uORFs 
(upstream  ORFs)  [16].  The  shorter  ER/32  3-UTR  (^2-short) 
overcame  this  inhibitory  influence,  whereas  the  longer  form 
(/32-long),  which  included  the  entire  sequence  of  /32-short, 
had  little  effect.  Regulatory  elements  in  both  5'-  and  3-UTRs 
have  important  roles  in  determining  translational  efficiencies, 
facilitated  by  transcript  circularization  induced  by  interaction  of 
poly(A)-binding  protein  binding  at  the  3'  end  of  transcripts  with 


eIF4G  (eukaryotic  initiation  factor  4G)  at  5'  ends  [33,34].  /32- 
short  apparently  interacted  with  5-UTRa,  resulting  in  a  loss  of 
inhibitory  structures  and/or  uORF  translation,  whereas  the  same 
sequences  were  unable  to  interact  in  this  way  in  the  context  of  /32- 
long.  Therefore  differential  use  of  polyadenylation  sites,  defining 
3-UTR  length,  exerts  translational  control.  It  is  interesting 
that  faster-proliferating  cells  and  cancer  cells  preferentially  use 
proximal  polyadenylation  sites  [35,36];  in  the  case  of  ER/32,  this 
would  lead  to  increased  translation  that  may  explain  how  ER/32 
is  up-regulated  in  breast  cancer  cells  [10,14]  despite  promoter 
methylation  and  down-regulation  of  ER/3  transcripts  [29,37]. 
Shortening  of  3-UTRs  has  been  associated  with  increased  protein 
expression  for  a  number  of  cancer-related  genes  [36].  For  these 
genes,  this  resulted  from  increased  mRNA  stability  and/or  loss 
of  microRNA-mediated  translational  repression  when  compared 
with  their  longer  3-UTRs.  For  ER/32,  this  is  clearly  not  the  case, 
since  ^2-short  dominantly  derepressed  the  influence  of  5-UTRa, 
as  opposed  to  merely  lacking  repressive  3'  elements  present  in 
/32-long.  ER/3 5  3-UTRs  also  tended  to  specify  relatively  efficient 
translation,  but,  in  this  case,  only  when  modifying  the  influences 
of  5 -UTRc  or  UTR-E1.  In  MCF7  cells,  this  was  only  apparent  for 
the  ER/3 5  3-UTR  that  was  actually  cloned  from  these  cells.  As  for 
ER/32,  the  data  do  not  support  a  simple  model  of  shorter  3-UTRs 
having  fewer  repressive  elements,  since  both  lengths  of  ER/3 5 
3-UTR  acted  to  derepress  the  influences  of  inhibitory  5-UTRs. 

Finally,  we  investigated  influences  of  the  ER  ligand  E2 
on  transcriptional  activity  of  ER/3  promoters  and  translation 
specified  by  ER/3  UTRs.  We  found  E2  exerted  relatively  mild, 
although  statistically  significant,  influences  at  multiple  levels, 
with  potential  to  change  the  balance  of  ER/3  isoform  expression. 
Unfortunately,  we  were  unable  to  confirm  whether  this  balance  is 
altered,  since  quantitative  detection  of  endogenous  ER/3  protein 
isoforms  in  cell  lines  is  unreliable  with  the  antibodies  currently 
available.  First,  we  found  E2  induced  transcription  in  MCF7  cells 
only  with  the  OK  promoter  (Figure  6),  a  promoter  we  have  shown 
previously  to  exhibit  a  relative  preference  for  producing  ER/31 
and  ER/3  5  transcripts  in  this  cell  type  [16].  Secondly,  we  found 
E2  induced  a  wide  range  of  changes  in  translational  efficiencies 
of  transcripts  with  ER/3  UTR  pairs  (Figure  7);  the  overall  effect 
in  any  cell  type  would  be  defined  by  the  relative  proportions  of 
each  UTR  pairing  within  the  total  pool  of  transcripts  for  that 
isoform.  E2  is  known  to  modulate  levels  of  microRNAs  [38,39], 
thus  these  constitute  potential  mediators  of  its  influence  on  ER/3 
translational  efficiencies.  Indeed,  we  have  found  that  miR-92  can 
regulate  expression  of  ER/31  acting  at  its  3-UTR  [40].  However, 
we  find  that  different  UTR  pairs  respond  differently,  therefore 
we  infer  that  accessibility  of  individual  microRNA-binding  sites 
in  the  context  of  interactions  between  5'-  and  3-UTRs  would 
influence  whether  particular  transcripts  respond  to  changes  in 
levels  of  any  particular  microRNA. 

In  conclusion,  we  have  revealed  novel  mechanisms  controlling 
ER/3  expression  that  help  explain  the  reported  lack  of  concordance 
between  ER/3  mRNA  and  protein  levels,  and  the  differential 
expression  of  ER/31,  ER/32  and  ER/35,  and  give  new  insights 
into  the  regulation  of  ER  function.  To  the  best  of  our  knowledge, 
this  is  the  first  time  cross-talk  between  multiple  5'-  and  3-UTRs 
has  been  implicated  in  the  differential  regulation  of  translation  of 
different  protein  isoforms  from  one  gene,  demonstrating  how  the 
ER/3  gene  provides  a  model  for  study  of  complex  gene  regulatory 
pathways. 
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Table  SI  Primer  sequences 


Primer  Sequence  (5'-*3') 


GFP  reporter  construction 
5'-UTRE1  + 

5'-UTR  - 
ER/31  3'-UTR  + 

ER/31  3'-UTR  - 
ER>02  3;-UTR  + 

ER>02  3'-UTR  short  - 
ER>02  3'-UTR  long  — 

ER/J5  3'-UTR  + 

ER£5  3'-UTR  short  - 
ER>05  3'-UTR  long  — 

Luciferase  reporter  construction 
Promoter  OK  + 

Promoter  OK  - 
Promoter  ON  + 

Promoter  ON  - 
Promoter  El  (5'  section)  + 
Promoter  El  (5'  section)  - 
Promoter  El  (3'  section)  + 
Promoter  El  (3'  section)  - 
5'-RACE 
First  round 
Second  round 
3'-RACE 
ER/31  first  round 
ER/31  second  round 
ER>02  first  round 
ER>02  second  round 
ER>05  first  round 
ER>05  second  round 
ER/J5  third  round 
Modified  adapter  primer  for  ER^2 
Real-time  PCR 
UTRa  + 

UTRc  + 

UTRa/c  - 
UTRE1  + 

UTRE1  - 
GFP  + 

GFP  — 

36B4  + 

36B4- 
Total  ER/3  + 

Total  ER/3  - 
ER/31  + 

ER/31  - 
ER>02  + 

ER>02  - 
ER>05  + 

ER>S5  — 


AGCAGAGCTCTCTGGCTAACTAGAGAACCCACCATTATACTTGCCCACGAATC 

AGCACTCGAGGTCTTGAGATAACAGCTGAGAAAACACC 

ATGCGGATCCCGCCTGGCCCTGAGG 

ATGCAAGCTTGAGGCCATTGAGTGTGGAAAC 

ATGCGGATCCCCCTCTAATCAACTCGG 

ATGCAAGCTTAAAGTATTTTAACTCTTTCTTTAAAA 

ATGCAAGCTTGCCCATTTAAGTCCAGTAGC 

ATGCGGATCCGGAGCTGCTCTGCTTG 

ATGCAAGCTTACACTTTTCCCAAATCAC 

ATG  C  AAG  CTTTTTTTTTT  C  ATG  G  ATTAC  AAT 

TATGCGGTACCGGCAATCAGATCCATTAGTTAAGC 

TATGCGCTAGCGTGGACGCCTACGAGGAGGGAGCG 

TATGCGGTACCTCCTTTCTCACTAGGTGGCTTCAGG 

TATGCGCTAGCGGGTGTCCAAAAAGCCAGCAGC 

TATGCGGTACCTCAAGACGAGCCTGGCCAACATGG 

TATGCGCTAGCCATTATAATGTTCTCAAAGATTCG 

ATGCGGTACCTAGAGTTTGGTGGGAAGCAGAGAGCTG 

ATGCGGTACCTCCCAAAGTGCTGGGATTACAGGCG 

CGTGCTCCAGGGGTAAGAT 

CAGGGGTAAGATGGATTGACTGC 

TGGAACATCTGCTCAACATGAAGTGC 

AATGTGGTCCCAGTGTATGACCTGC 

GGCTAACCTCCTGATGCTCCTGTCC 

GTCAGGCATGCGAGGGCAG 

GGCTAACCTCCTGATGCTCCTGTCC 

GTCAGGCATGCCAGGTACGC 

CGCCCTAAGGAGCTGCTCTGC 

GGCCACGCGTCGACTAGTACTTTTTTTTTTTTTTTTT 

AGTTACTGAGTCCGATGAATGTGCTTG 

CGGGAGACCCCCCCTAATGC 

CTCAAAGATTCGTGGGCAAGTATAATG 

CTTGGCAATATCTTCTGTGTTTCTTTACAG 

GGCTAGATGGTGAGTTTTTTATATCCATGT 

CAACCACTACCTGAGCACCCAGTC 

GGCGGCGGTCACGAACTC 

GAAACTCTGCATTCTCGCTTCC 

GATGCAACAGTTGGGTAGCCA 

TGGGCACCTTTCTCCTTTAGTGG 

GCTTCACACCAGGGACTCTTTTGAG 

CGCCTGGCTAACCTCCTGATG 

GAGCAGATGTTCCATGCCCTTG 

CGCCGTGACCGATGCTTTG 

CCTTTTCTGCCCTCGCA 

GCATCTCCTCCCAGCAGCAATC 

GCACATAATCCCATCCCAAGCC 
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